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Abstract 

Traditionally  flight  control  systems  have  used  linearized  equations  of  motion 
solved  around  a  single  trim  point.  This  thesis  proposes  a  nested- loop  controller  directly 
solved  from  the  equations  of  motion. 

The  control  equations  were  developed  as  a  solution  to  asymmetrically  trimmed 
flight  conditions.  A  two- loop  design  was  proposed  for  the  controller.  The  outer  loop 
modeled  the  aircraft  as  a  point  mass  and  all  forces  were  balanced  to  find  the  aircraft 
states.  The  equations  input  the  control  variables  and  output  the  aircraft  states.  The  inner- 
loop  utilizes  the  six- degree  of  freedom  model  of  the  aircraft  to  solve  the  moment 
equations.  With  the  input  states,  the  required  control  surface  deflections  are  calculated. 

The  control  equations  were  investigated  for  typical  flight  conditions  to  find  the 
predicted  aircraft  control  settings.  The  control  equations  were  implemented  using 
aeromodel  data  for  aEearjet-25.  The  aeromodel  data  was  updated  in  flight  test.  The 
predictions  from  the  control  equations  were  then  compared  to  flight  test  results.  The 
model  was  able  to  predict  the  required  elevator  deflection  for  simple  longitudinal  cases  in 
level  and  climbing  flight  to  within  tolerances.  The  simple  lateral- directional  cases  were 
not  as  accurate  as  the  longitudinal  investigations.  As  complex  maneuvers  were 
investigated,  the  model  predictions  did  not  match  the  flight  test  results.  The  complex 
maneuvers  were  not  reproduced  in  flight  test  to  match  the  flight  parameters  calculated 
with  the  model.  Also  the  lateral- directional  stabihty  derivatives  and  measurements  had 
larger  errors  than  the  longitudinal  variables. 
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OUTER- LOOP  CONTROL  IN  ASYMMETRICAL  TRIMMED  LLIGHT  CONDITIONS 


I.  Introduction 


Background 

Automatic  flight  control  systems  are  a  vital  part  of  most  aircraft  designs.  The 
autopilot  controls  both  the  longitudinal  and  lateral- directional  motion  of  the  vehicle  to 
perform  the  desired  mission.  The  control  system  decreases  pilot  workload,  enhances 
flight  safety,  and  increasingly  is  used  to  .stabilize  aircraft  that  otherwise  would  not  be 
stable.  Also,  with  the  increased  use  of  Unmanned  Aerial  Vehicles  (UAVs),  the  automatic 
flight  control  system  is  vital  to  achieve  the  required  performance.  Traditionally, 
automatic  flight  control  systems  have  used  hnearized  equations  of  motion  to  simplify 
aircraft  motion  around  an  equihbrium  point.  Linearization  methods  have  been  studied 
extensively  and  are  implemented  on  most  modem  aircraft  (1:  Chapter  5).  It  is  proposed 
to  solve  the  aircraft  equations  of  motions  to  find  the  necessary  control  surface  deflections 
to  achieve  a  desired  flight  path  and  orientation.  The  developed  equations  are  applicable 
to  most  of  the  flight  envelope  and  not  fixed  to  the  equihbrium  point. 

Lrom  the  time  of  the  Wright  brothers,  aircraft  stability  and  control  has  been  a 
central  concern  of  aircraft  design.  To  meet  design  performance  requirements,  a  balance 
must  be  achieved  between  aircraft  stabihty  and  maneuverabihty.  Increasingly,  modem 
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aircraft  are  designed  to  be  inherently  unstable  to  meet  increased  operational  demands. 
The  F- 16  fighter  aircraft  is  unstable  in  the  unaugmented  airframe  configuration,  and 
requires  an  automatic  flight  control  system  to  maintain  steady  flight  (2:  Class  2). 

Automatic  flight  control  systems  have  been  studied  extensively,  and  many  schemes  have 
been  designed  (3:  367,  386). 

The  traditional  design  method  for  flight  control  systems  is  to  assume  flight  at  a 
stable  equihbrium  point  with  small  perturbations  from  that  point.  The  aircraft  equations 
of  motion  are  linearized,  and  small  angle  assumptions  are  normally  made  to  simphfy  the 
equations.  This  method  is  generally  only  vahd  for  a  small  region  and  must  be 
recalculated  if  flight  conditions  change  from  the  equihbrium  point.  This  requires  the 
automatic  flight  control  system  to  be  divided  into  smah  regions  of  apphcabihty  (gain 
scheduhng).  Changing  from  one  gain  region  to  another  may  cause  discontinuities  in  the 
equihbrium  point  and  jumps  in  the  commanded  control  inputs.  Therefore,  it  is  desirable 
to  design  an  automatic  flight  control  system  that  is  apphcable  for  most  of  the  flight 
envelope. 

Overview 

A  design  is  envisioned  to  solve  the  required  control  surface  deflection  angles  from 
the  aircraft  equations  of  motion  for  a  desired  flight  condition.  A  design  was  proposed 
utilizing  a  nested  loop  confroher  to  solve  the  equations  of  motion  (reference  4).  The 
outer- loop  dynamics  of  the  aircraft  were  modeled  as  a  point  mass,  using  nonhnear 
equations  of  motion.  The  flight  state  variables  (velocity  (V),  ghde  path  angle  (y)  and 
course  angle  (H)  change  rate)  were  nondimensionalized  and  used  to  describe  the  desired 
trimmed  steady  state  flight  trajectory.  The  flight  trajectory  is  defined  as  the  path  the 
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aircraft  center  of  gravity  travels  through  the  atmosphere.  Additionally,  to  solve  the 
equations  either  sideshp  angle  ((3)  or  bank  angle  ((|)v)  must  be  specified  for  the  given 
flight  condition.  Equations  were  developed  to  output  angle  of  attack  (a),  thmst  setting 
(T),  sideslip  angle,  and  bank  angle  from  the  state  variable  inputs. 

The  remainder  of  the  design  was  developed  and  implemented  for  verification  on 
an  actual  aircraft.  The  outer- loop  outputs  were  dimen sionalized  for  the  aircraft  model 
data.  The  results  of  the  outer- loop  equations  were  in  the  wind  axis,  which  were  then 
transformed  to  the  body  axis  to  become  the  inputs  into  the  inner- loop  equations.  The 
inner- loop  equations  use  the  calculated  a,  T,  (3,  (|)v  and  specified  flight  orientation  to  find 
the  required  aileron  (5a),  elevator  (5e),  and  rudder  (5r)  surface  deflection  angles.  The 
flight  orientation  is  defined  as  the  aircraft  rotation  and  angular  acceleration  about  the 
center  of  gravity.  As  long  as  the  required  deflections  are  smaller  than  the  deflection 
limits,  the  desired  trimmed  flight  path  can  be  achieved  by  commanding  the  control 
surfaces.  The  control  laws  were  implemented  using  a  Leaijet-25  model  to  determine  the 
validity  of  the  approach.  Finally,  the  trim  condition  parameters  were  estabhshed  in  flight 
tests  to  verify  the  predicted  results. 

The  remainder  of  the  paper  is  in  the  following  format.  In  Chapter  2,  the  three 
axes  used  for  the  equations  of  motions  (wind  axis,  navigation  axis,  and  body  axis)  are 
introduced.  The  Learjet-25  aircraft  model  is  also  presented  and  all  model  parameters  are 
hsted  in  Table  1.  Next,  the  equations  previously  derived  to  calculate  the  outer- loop  are 
reviewed  and  modified  to  match  the  defined  axes.  The  aircraft  trim  equations  for  the 
inner- loop  dynamics  are  derived  in  Chapter  3.  In  Chapter  4,  selected  flight  trajectories 
are  investigated  for  trimmed  equUibrium  flight  conditions.  First,  straight  and  level  flight 
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is  explored.  A  wings  level  climb  and  wings  level  descent  are  then  calculated.  Next, 
crabbing  flight  (sideshp  angle  not  equal  to  zero)  is  determined.  Steady,  level  turn  and 
corkscrew  turn  (cUmbing  and  descending)  maneuvers  are  then  investigated.  Next,  a 
barrel  roll  maneuver  (constantly  changing  roll  angle)  is  reproduced.  Then  the  AC  - 130 
gunship  level  orbit,  firing  maneuver  is  examined.  Finally,  a  steady  pitch- over  maneuver 
is  modeled.  In  Chapter  5,  the  effect  of  errors  in  the  stabihty  derivatives  are  investigated 
by  varying  the  known  derivative  values  ±40%  to  determine  the  change  in  the  control 
surface  predictions.  Control  surface  deflection  changes  greater  than  1.0  degree  are 
considered  to  have  a  significant  impact  on  the  equation  predictions.  In  Chapter  6  the 
predictions  from  Chapter  4  were  vahdated  by  flight  test  data.  Selected  trim  conditions 
were  examined  in  the  flight  test  program  and  compared  to  the  theoretically  calculated 
values.  Chapter  7  includes  the  conclusions  and  recommendations. 
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II.  Model  Definition 


Overview 

In  Chapter  2,  the  point  mass  model  is  developed  to  provide  a  basis  for  the  control 
laws  for  the  trajectory  of  the  aircraft.  The  outer- loop  equations  are  found  from  the 
nonhnear  equations  of  motion  and  apphed  to  the  actual  Leatjet-25  parameters.  The 
outer- loop  model  equations  are  reviewed  to  provide  a  basis  for  the  rest  of  the  study.  The 
equations  were  previously  developed  (reference  4),  but  are  reviewed  for  clarity  due  to  the 
need  to  redefine  the  reference  axes.  The  control  equations  are  developed  from  balancing 
the  forces  and  moments  (trimmed  flight).  In  Chapter  3,  the  outer- loop  equation  outputs 
are  used  with  the  transformation  matrices  and  the  remainder  of  the  equations  of  motion  to 
find  the  control  surface  deflections  required  to  achieve  the  desired  flight  trajectory  and 
orientation.  In  Chapter  4,  the  equations  developed  in  the  previous  two  chapters  are  used 
to  simulate  selected  flight  trajectories.  By  choosing  the  correct  input  states,  the  required 
control  surface  deflections  for  the  desired  flight  path  are  found.  The  impacts  of 
variations  in  the  stability  derivatives  are  investigated  in  Chapter  5  to  determine  which 
parameters  require  the  most  accuracy.  Finally,  in  Chapter  6  the  theoretical  values  tested 
using  the  Learjet-25  in  flight  tests  are  compared  to  the  theoretical  calculations. 

Axes  Definition 

Background  information  and  the  aircraft  model  are  now  presented.  First,  the 
pertinent  reference  frames  are  introduced,  including  the  coordinate  transformations.  The 


5 


outer- loop  equations  are  re- derived  including  the  necessary  changes  for  the  redefined 
axes.  Finally,  the  Leaijet-25  aerodynamic  parameters  are  provided. 

Three  reference  frames  are  used  in  the  development  of  the  control  equations:  body 
axes,  wind  axes  and  the  navigation  frame.  The  body  axes  are  affixed  to  and  move  with 
the  aircraft.  The  body  axes  parameters  are  easiest  to  measure  since  they  remain  constant 
with  respect  to  the  aircraft.  The  wind  axes  are  fixed  to  the  aircraft  center  of  gravity  (eg) 
and  utilize  the  free- stream  airflow  as  a  reference.  Wind  axes  reference  frames  are  used 
extensively  in  analyzing  flight  mechanics.  The  navigation  frame  is  attached  to  the  earth 
and  provides  an  approximate  inertial  reference  frame.  The  axes  definition  and 
relationships  between  reference  frames  follow. 

Body  Axes. 

The  body  axes  frame  is  specified  in  relation  to  the  sfructure  of  the  aircraft.  The 
body  x-axis  (xb)  points  out  of  the  nose  of  the  aircraft.  The  body  y-axis  (yb)  is  orthogonal 
to  the  Xb-axis,  passing  out  the  right  wing,  and  the  body  z-axis  (zb)  completes  the  right- 
handed  coordinate  system,  pointing  in  a  downward  direction.  The  body  axes  system  is 
defined  relative  to  the  North- East- Down  (NED)  navigation  axes  frame  by  specifying  the 
3-2-1  Euler  angles:  Eirst,  the  heading  angle  (t|/)  is  the  rotation  about  the  navigation  z-axis 
(zn)  to  estabhsh  the  intermediate  1-axes.  Next,  the  pitch  angle  (0)  is  the  rotation  about 
the  yi-axis  to  find  the  intermediate  2- axes.  EinaUy,  the  roll  angle  ((|))  is  the  rotation  about 
the  X2-axis  to  produce  the  body  axes.  The  rotations  and  angles  are  shown  in  Eigure  1. 

The  inner- loop  calculations  are  performed  in  the  navigation  axes  frame. 
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Figure  1.  Navigation- to- Body  Axes  Transformation 


Navigation  Axes. 

The  navigation  axes  provide  an  approximate  inertial  reference  frame  and  are  used 
for  most  performance  calculations.  In  flight  mechanics  though,  the  wind  axes  frame  is 
typically  used.  The  navigation  axes  are  also  related  to  the  wind  axes  by  a  set  of  3-2- 1 
Euler  angles.  First,  the  course  angle  (H)  is  defined  as  the  rotation  about  the  ^-axis  to 
find  the  intermediate  3- axes.  Next,  the  flight  path  angle  (y)  is  defined  as  the  rotation 
about  the  ya-axes  to  find  the  intermediate  4- axes.  Finally,  (|)v  is  defined  as  the  rotation 
about  the  X4-axis  to  produce  the  wind  axes.  The  Euler  angles  H,  y,  and  (|)v  are  shown  in 
Figure  2. 
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Figure  2.  Navigation-to-Wind  Axes  Transformation 
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Wind  Axes. 


The  outer- loop  control  model  equations  use  the  wind- axes  system.  The  wind  axes 
are  defined  as  the  Xw-axis  pointing  along  the  forward  velocity  vector  of  the  free- stream 
airflow.  The  Zw-axis  is  defined  negative  along  the  aircraft  hft  vector,  perpendicular  to  the 
Xw-axis.  The  yw-axis  is  orthogonal  to  both  the  Xw-  and  Zw-axis,  pointing  to  the  right.  The 
relationship  between  the  wind  axes  and  the  body  axes  is  found  by  first  rotating  the 
aircraft  about  the  z^^axis  by  the  sideshp  angle  to  find  the  intermediate  5-axes.  The  wind 
axes  are  then  defined  as  the  rotation  of  the  angle  of  attack  about  the  new  ys-axis.  The 
rotations  and  angles  are  shown  in  Figure  3. 
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Figure  3.  Wind-to-  Body  Axes  Transformation 


Outer-Loop  Model  Development 

This  is  the  definition  of  the  aerodynamic  angles  used  in  this  paper.  Unfortunately, 
there  is  no  standard  definition  of  the  aerodynamic  angles,  a  and  (3  (1:62),  (5:36  and  88). 
Sideshp  angle,  (3,  is  normaUy  assumed  to  be  smaU,  so  this  is  not  normaUy  a  significant 
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issue.  For  large  angles  the  correct  definition  of  a  and  (3  is  important  and  a  rotation  order 
must  be  specified.  Sideslip  is  chosen  as  the  first  rotation  angle  because  in  a  wind-tunnel, 
the  model  would  normally  be  first  rotated  by  |3,  and  then  a  for  each  test  condition.  When 
the  aerodynamic  angles  are  small,  roll  angle,  (|)v,  is  approximately  equal  to  (|).  The  pitch 
angle,  0,  is  approximately  equal  to  the  sum  of  the  angles  a  and  y.  The  yaw  angle,  \\t,  is 
approximately  equal  to  the  combined  angle  of  H  and  (3 .  The  previous  outer- loop  study 
(reference  4)  had  different  definitions  for  the  axes,  but  with  the  above  definitions  of  the 
axes,  existing  aircraft  performance  data  could  be  used  without  modifications.  The 
previous  aircraft  equations  of  motion  were  altered  to  reflect  the  new  axes  definitions. 

Background. 

The  outer- loop  equations  were  previously  derived  for  a  generic  fighter  aircraft. 
The  derivation  of  the  equations  of  motion  is  reviewed  with  the  new  definition  of  the 
coordinate  axes.  Even  though  the  yr  and  Zb-  axes  were  rotated  180  degrees,  the  only 
difference  in  the  equations  is  the  definition  of  the  fiight  path  angle,  y,  which  is  now 
negative  with  the  new  definitions.  The  side- force  (Fy)  has  the  same  sign,  but  now  is 
defined  as  positive  to  the  right.  The  forces  acting  on  the  aircraft  are  shown  on  Figure  4. 
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Point-Mass  Model  Definition. 


y-i 


Figure  4.  Force  Definitions 


The  aerodynamic  forces  are 


2  ^La 

Vo 

(1) 

where  q  =  ^nV^^S^  (dynamic  pressure) 

D=q^(CD,a+Ci,J 

(2) 

-  S 

Fy=q  —  ^  P 

y  4  2  s  ya 

(3) 

T=mgp  (thrust) 

(4) 

The  non-dimensionalized  thmst  setting  is  p. 
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The  aircraft  equations  of  motion  using  wind  axes  are 


The  same  assumptions  and  simphfications  of  the  previous  study  (reference  4)  to 
the  equations  were  applied  to  find  the  control  variables  |i,  a,  |3,  and  (|)v  as  functions  of  the 
state  variables  (V,  y  and  H).  Details  of  the  calculations  are  found  in  Appendix  A.  With 
the  equations  of  motion  derived,  specific  trim  conditions  are  then  examined.  By 
describing  the  trimmed  flight  condition  with  respect  to  the  state  variables,  control  settings 
were  obtained  for  the  trim  conditions. 

The  stabihty  derivative  was  used  for  the  Learjet-25  model.  Often  the  drag  is 

2  2 

defined  as  a  function  of  (Cl^  a  ),  but  the  available  model  used  .  The  model 
will  not  be  as  accurate  as  the  trim  condition  gets  farther  from  flight  condition  where 

was  determined  since  does  not  account  for  the  second-order  effects  of  Cl^  ^ .  Using 

2 

may  increase  the  accuracy  away  from  the  point  where  the  stabihty  derivatives  were 
calculated. 
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Straight  and  Level  Flight. 


For  trimmed  flight,  the  derivatives  of  the  state  variables  5<,  and  1^)  are  set 
equal  to  zero.  Details  of  the  calculations  are  found  in  Appendix  A.  Again,  one  of  the 
control  variables  (3  or  (|)v  must  be  chosen  for  the  specified  trim  condition. 

The  trim  control  law  for  a  specified  trim  (3  is 


sin(|)v  = 


Cy,S,q 

X"nmgS„  cos  a 


(8) 


COS(|)v  = 


1 


1- 


^  C,,S.qV^a 
V„^nmgS„  cos  a 


(9) 


a- 


C, 


X  2  A 

Vjj  nmg  COS  a 

qv^ 


-a. 


(10) 


^‘=-^(CD.-C„.a^)+. 

mgV„ 


f Cq  ncosa^ 

2 

Ca.  J 

+  nsin  a 


(11) 


and  the  trim  control  law  for  a  specified  trim  (|)v  is 


p^V^nmgS, 

C.S.q  V=  ■ 


a  = 


1  (  V^^nmg  cos  a 


C, 


c 

D„  / 

u  = - (n 

C 

^4. 


qV^ 


cos  acos( 


cosd^ 


qV^ 


mgV„ 


(12) 

(13) 

(Cd^  )  +  nsin  a 

(14) 
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Steady  Turn  Flight. 


For  steady  turning  flight,  ’^and  5<  are  set  equal  to  zero,  but  1^  is  a  constant 
specifled  course  change  rate,  05,j, .  Therefore,  H(t)  is  equal  to  t.  For  a  specified  trim 
P,  the  control  law  is 


(18) 


Steady  Rolling  Flight. 

For  steady  roUing  flight,  ^  and  S<  are  set  equal  to  zero,  but  the  velocity 
vector  roll  rate  is  a  constant  specifled  rate,  05^  .  Therefore,  (|)v(t)  is  equal  to  05^  t. 
Again,  for  a  specifled  05^ ,  the  trim  control  law  is 
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From  equation  (12) 


|3(t)= 


VotimgS 

Cy.S.q 


cos  a  .  _ 


(19) 


From  equation  (13) 


a(t)= 


C, 


V^nmg  cosacosco^t 


(20) 


From  equation  (14) 


F(t> 


— ^ncos  acosCGLt  +  nsin  a 


(21) 


Learjet-25  Aeromodel. 

The  aircraft  model  used  is  the  Learjet-25  aircraft  operated  by  Veridian  for  United 
States  Air  Force  Test  Pilot  School  (TPS).  The  aircraft  data  is  listed  in  Table  1  (reference 
7). 
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Table  1.  Learjet  Aircraft  Data 


AIRCRAFT  PARAMETER 


Weight  (Standard  Fuel  Load) 


Wing  Area 


Wing  Span 


Mean  Aerodynamic  Chord 


Vertical  Tail  Area 


Aircraft  Lift-Curve  Slope 


Aircraft  Zero- Lift  Coefficient 


Lift  Coefficient  due  to  Elevator 


Lift  Coefficient  due  to  Pitch  Rate 


Aircraft  Drag  Coefficient  due  to  Lift 


Aircraft  Zero- Lift  Drag  Coefficient 


Drag  Coefficient  due  to  Pitch  Rate 


Mass  Moment  of  Inertia  (xb-axis) 


Mass  Moment  of  Inertia  (yb-axis) 


Mass  Moment  of  Inertia  (zb-axis) 


Product  of  Inertia  (xb-  and  Zb-axes) 


Side- Force  Coefficient  due  to  Sideshp 


Side- Force  Coefficient  due  to  Rudder 


Side- Force  Coefficient  due  to  Aileron 


Side- Force  Coeff.  due  to  Roll  Rate 


Side- Force  Coeff.  due  to  Yaw  Rate 


Roll  Moment  Coeff.  due  to  Sideshp 


Roll  Moment  Coeff.  due  to  Rudder 


RoU  Moment  Coeff.  due  to  Aileron 


RoU  Moment  Coeff.  due  to  RoU  Rate 


Roll  Moment  Coeff.  due  to  Yaw  Rate 


Yaw  Moment  Coeff.  due  to  Sideshp 


Yaw  Moment  Coeff.  due  to  Aileron 


Yaw  Moment  Coeff.  due  to  Rudder 


Yaw  Moment  Coeff.  due  to  Roh  Rate 


Yaw  Moment  Coeff.  due  to  Yaw  Rate 


VARIABLE 


W 


VALUE 


13,500 


231.8 


39.5 
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76.5 


4.967 


0.0105 


0.3631 


37.357 


0.1146 


0.0235 


4.927 


20,000 


22,900 


40,000 


1,980 


-0.7620 


0.1423 


-0.0613 


8.549 


47.423 


-0.0945 


0.0198 


-0.0791 


-25.573 


22.221 


0.0830 


-0.0011 


-0.0516 


-3.676 


-6.666 


UNITS 


lbs 
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Table  1.  (cont.) 


AIRCRAFT  PARAMETER 

VARIABLE 

VALUE 

UNITS 

Pitch  Coeff.  due  to  Angle  of  Attack 

-0.8921 

rad' 

Pitch  Coefficient  due  Alpha  Change 

C.. 

-5.615 

s/rad 

Zero- Lift  Pitch  Coefficient 

0.0529 

Pitch  Coefficient  due  to  Pitch  Rate 

-864.347 

raffi^ 

Pitch  Coefficient  due  to  Elevator 

-0.8434 

raffi^ 

Maximum  Thmst  Right  Engine 

Ti 

^max 

5000 

lbs 

Maximum  Thmst  Eeft  Engine 

T 

5000 

lbs 

X- direction  Distance  c.g.  to  Thmst 

XT 

0.1 

ft 

Y-direction  Distance  c.g.  to  Thmst 

yx 

3 

ft 

Z-direction  Distance  c.g.  to  Thmst 

Zt 

0.1 

ft 

Engine  Thmst  Angle 

^  j 

0.0 

rad 

Maximum  Elevator  Deflection  Angle 

^max 

0.2618 

rad 

Minimum  Elevator  Deflection  Angle 

^min 

-0.2618 

rad 

Maximum  Aileron  Deflection  Angle 

“max 

0.5236 

rad 

Minimum  Aleron  Deflection  Angle 

“mm 

-0.5236 

rad 

Maximum  Rudder  Deflection  Angle 

^max 

0.6109 

rad 

Minimum  Rudder  Deflection  Angle 

^min 

-0.6109 

rad 

Sea- level  Air  Density 

Po 

wkwiifagiriii 

Universal  Gas  Constant 

R 

1718 

ftlbs/slug°R 

Sea- level  Temperature 

To 

518.69 

“R 

Gravitational  Acceleration 

g 

32.174 

ft/sec^ 

Suimnary 

In  this  chapter  the  aircraft  model  was  first  described.  New  axes  definitions  were 
presented  and  the  relationships  between  body,  wind  and  navigation  axes  were  detailed. 
The  outer- loop  equations  of  motion  were  then  redefined,  to  produce  new  control  laws. 
With  the  new  control  laws,  state  variable  values  were  chosen  to  describe  a  desired  flight 
trajectory,  and  the  resulting  control  settings  were  determined.  Therefore,  the  outer  loop 
model  uses  the  state  variables  (V,  y,  and  H)  and  a  specified  value  for  either  |3  or  (|)v,  and 
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outputs  the  necessary  control  variables  (a,  |3,  (|)v,  and  |i).  The  inner- loop  equations  use 
the  aircraft  state  values  and  the  above  control  settings  to  determine  the  aircraft  control 
surface  deflections  (5a,  6e,  and  5r)  in  the  next  chapter.  Finally,  the  Learjet  aircraft 
parameters  necessary  to  calculate  the  inner-  and  outer- loops  were  provided. 
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III.  Inner-Loop  Analysis 


Overview 

In  Chapter  2,  the  point  mass  model  was  developed  to  provide  a  basis  for  the 
control  laws  for  the  trajectory  of  the  aircraft.  The  outer- loop  equations  were  found  from 
the  nonlinear  equations  of  motion  and  apphed  to  the  actual  Learjet-25  parameters.  In 
Chapter  3,  the  outer- loop  equation  outputs  are  used  with  the  transformation  matrices  and 
the  remainder  of  the  equations  of  motion  to  find  the  control  surface  deflections  required 
to  achieve  the  desired  flight  trajectory  and  orientation.  With  the  outer- loop  equations 
developed  in  the  last  chapter,  the  inner- loop  equations  are  established  to  provide  aircraft 
control  laws  for  the  Learjet-25.  Selected  trim  conditions  are  then  examined  to  determine 
the  equations  to  predict  the  required  control  surface  deflections.  In  Chapter  4,  the 
equations  developed  in  the  previous  two  chapters  are  used  to  simulate  selected  flight 
trajectories.  By  choosing  the  correct  input  states,  the  required  control  surface  deflections 
for  the  desired  flight  path  are  found.  The  impacts  of  variations  in  the  stabihty  derivatives 
are  investigated  in  Chapter  5  to  determine  which  parameters  require  the  most  accuracy. 
Finally,  in  Chapter  6  the  theoretical  values  tested  using  the  Learjet-25  in  flight  tests  are 
compared  to  the  theoretical  calculations. 

Inner-Loop  Model  Development 

The  inner- loop  model  uses  the  trim  state  variables  and  the  trim  control  variables 
from  Chapter  2  to  determine  the  needed  trimmed  aircraft  control  surface  deflections. 
First,  the  navigation  axes  Euler  angles  are  calculated  from  the  known  wind  axes  using  the 


18 


relevant  transformation  matrix.  With  the  Euler  angles  and  the  specified  angular  rotations 
for  the  desired  flight  trajectory  and  orientation,  the  aircraft’s  angular  rates  resolved  in  the 
body  axes  are  determined.  Next,  with  the  six- degree- of- freedom  aircraft  equations  of 
motion,  a  steady- state  trim  condition  was  chosen  to  find  the  required  control  surface 
deflections  from  the  known  aircraft  parameters.  Finally,  assumptions  were  made  for  the 
aircraft  to  simphfy  the  control  surface  deflection  equations  for  the  specific  trim 
conditions. 

The  body,  wind  and  navigation  axes  were  defined  in  Chapter  2,  with  the  required 
rotation  angles  to  translate  between  coordinate  systems.  The  angles  (a,  (3,  (|)v,  Y,  and  H) 
are  used  in  the  outer- loop.  The  inner- loop  equations  use  the  Euler  angles  (t|/,  0,  and  (|)). 
The  relevant  axes  conversion  matrix  is  used  to  determine  the  Euler  angles. 

Navigation  Axes  Rotation. 

First,  the  navigation- to- body  axes  conversion  matrix  was  determined  using  Figure 
1.  The  following  three  rotations  define  the  axes  rotations 


cos'E  - 

-sin'E  0 

C0  = 

sin  'E  cos  'E  0 

0 

0  1 

cose 

0  sin  e 

C,= 

0 

1  0 

-sin  e 

0  cose 

'1  0 

0 

C,= 

0  costf) 

-sin  (|) 

0  sintf) 

COStf) 
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The  coordinate  transformation  matrix  is  given  by 

r"  =r  r  -r 

sin(|)  sin  ecos  'P  -  cos  (|)  sin  'P  cos  (|)  sin  ecos'P  +  sin  (|)  sin  'P 
sin(|)  sin  esin'T+ cos(|)  cos'T  cos  (|)  sin  esin 'T  -  sin  (|)  cos 'T  (18) 
sin  (|)  cos  e  cos(|)  cos  e 


Wind  Axes  Rotation. 

The  wind-to-body  axes  rotations  are  shown  in  Figure  2.  The  following  two 
transformations  define  the  axes  rotations 

cos  a  -  sin  a  0 
C.  =  sin  a  cosa  0 

a 

_  0  0  1 

cosa  0 
C,  =  0  1 

a 

sin  a  0 

where  the  transformation  matrix  is 

cos  acosa  -  sin  a 
C"  =  C.  •  =  cos  asina  cos  a 

-sin  a  0 


cos  ecos'F 
C|5  =  cosesin'T 
-sin  e 
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Wind  to  Navigation  Axes  Rotation. 


The  wind-to- navigation  axes  rotations  are  shown  in  Figure  3.  The  following  three 
transformations  define  the  axes  rotation 

cosH  -sin  H  0 
Cjj  =  sin  H  cos  H  0 
_  0  0  1 

cos  a  0  sin  a 
C,  =  0  10 

a 

sin  a  0  cos  a 

'10  0 
=  0  cos(|)^  -sin(|)^ 

_0  sintf)^  costf)^ 

where  the  transformation  matrix  is 
C"  =  C„  C-  •  Q 

w  a  <|)v 

cosacosH  sm(|)^sin  acosH-cos(|)^sin  H  cos(|)^  sin  acosH  +  sin  (|)^  sin  H 
=  cos  asinH  sin  (|)y  sin  asinH  +  cos(|)y  cos  H  cos(|)  sin  asinH  -  sin  (|)^  cos  H  (24) 
-sin  a  sin(|)^cosa  cos(|)yCosa 

Combined  Rotations. 

Now  the  combined  wind  and  body  axes  rotation  is  given  by 


(25) 
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cos  ecos'P 
cosesin'P 
-sin  e 


sin  (|)  sin  Ceos'?  -  cos  (|)  sin  *?  cos  (|)  sin  coos'?  +  sin  (|)  sin  '? 
sin(|)  sin  esin '?  +  cos(|)  cos '?  cos  (|)  sin  esin'?  -  sin (|)  cos '? 
sin  (|)  cos  e  cos(|)  cos  e 


cos  acosH  sin(|)^sin  acosH-cos(|)^,  sinH 
cos  asinH  sin(|)^  sin  asinH  +  cos (|)^  cos  H 
-sin  a  sin(|)^cosa 


cos  (|)^,  sin  acosH+  sin  (|)^  sin  H 
cos  (|)^  sin  asinH  -  sin  (|)^  cos  H 
cos(|)^cos  a 


cos  acosa 

-sin  a 

sin  acosa 

C12 

C13’ 

cos  asina 

cosa 

sin  asina 

= 

^21 

^22 

^23 

-sin  a 

0 

cosa 

.C31 

^32 

^33  _ 

where 


Cj  j  =  cos  a  cos  Hcos  a  cos  a  +  (sin  (f)^  sin  a  cos  H  -  cos  (f)^  sin  H)  cos  a  sin  a  + 
-  (cosef)^  sin  acos  H  +  sin  (j)^  sin  H)  sin  a 
Cj2  = -cos  acos  H  sin  a  +  (sin  (f)^  sin  a  cos  H -cos  (|)^  sin  H)cosa 
Cj3  =  cos  acos  H  sin  a  cosa  +(sin  (f)^  sin  a  cos  H  -cos  (j)^  sin  H)  sin  a  sin  a  + 

+  (cos  (j)^  sin  a  cos  H  +  sin  (f)^  sin  H)  cos  a 
C21  =  cos  asinH  cos  a  cosa  +  (sin  (f)^  sin  asinH  +  cosef)^  cos  H)  cos  a  sin  a  + 
-(cosef)^  sin  asin  H-sin  (l)^cosH)sin  a 
C22  =  -cos  asinH  sin  a  +  (sin  (|)^  sin  asinH  +  cos (j)^  cos  H)  cos  a 
C23  =  cos  asinH  sin  a  cos  a  +  (sin  (j)^  sin  asinH  +  coscj)^  cos  H)  sin  a  sin  a  + 

+  (cos(|)^  sin  asin  H  -  sin  cos  H)  cos  a 

C3J  =  -sin  a  cos  a  cos  a  +  sin  (1\  cos  acosa  sin  a  -  cos  (j)^  cos  a  sin  a 
C32  =  sin  a  sin  a  +  sin  (f)^  cos  acos  a 

C33  =  -sin  asin  a  cos  a  +  sin  (|)^  cos  asina  sin  a  +cos(l)^  cos  acos  a 


Navigation  Angle  Determination. 

From  the  above  matrix  equations,  the  navigation  axes  angles  were  determined. 
Equations  for  both  the  sine  and  cosine  of  the  angles  were  determined  to  resolve  the 


quadrant  of  the  angle.  For  an  angle  between  0  and  —  both  sine  and  cosine  values  are 
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positive.  For  positive  sine  values  and  negative  cosine  values,  the  angle  is  between  —  and 


3^  ^ 

—  .  In  the  third  quadrant  (0  to  -  — ),  sine  values  are  negative  and  cosine  values  are 


^  3^ 

positive.  For  an  angle  between  -  —  and  -  —  both  sine  and  cosine  values  are  negative. 


The  navigation  axes  angles  were  found  by  solving  the  following  equations.  Setting 
individual  elements  of  the  transformation  matrix  equal 
The  first  row,  first  column  terms  of  both  matrices  are 

cosecos'F  =  Cjj 

The  second  row,  first  column  terms  of  both  matrices  are 

cos  e  sin  'F  =  C21 

Squaring  both  equations  and  adding  them  together 

cos^ecos^'F  +  cos^esin^'F  =  Cjj^  +  C2i^ 
which  with  cos^  'F  +  sin  ^  'F  =  1  yields 


cose 


e  =  Vc 


+  c 
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(26) 


From  the  third  row,  the  first  column  terms  are 


sin  e  =  -Cjj  (27) 
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From  the  first  row,  first  column  terms 


cos  'F  = 

cose 


(28) 


From  the  second  row,  first  column  terms 


sin'F  =  -^ 
cose 


(29) 


From  the  third  row,  third  column  terms 


costf)  = 


•"33 


cose 


(30) 


From  the  third  row,  second  column  terms 


sin  (f) : 


"32 


cose 


2  ,  2 


(31) 


Trim  Condition  Definitions 

Therefore,  two  expressions  were  found  for  each  of  the  navigation  axes  angles,  \|/, 
0,  and  (|),  as  functions  of  the  wind  axes  variables.  Specific  trim  conditions  are  then 
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examined,  to  determine  simplified  control  laws.  The  equations  are  not  defined  at  0  =  90°, 
therefore  the  flight  with  the  nose  straight  up  is  not  allowed. 

Trim  condition  I  corresponds  to  steady,  straight  and  level  flight.  A  straight  and 
steady  climb  is  described  by  trim  condition  2.  Trim  condition  3,  a  straight  and  steady 
descent,  is  the  opposite  of  trim  condition  2.  In  trim  condition  4,  level  crabbing  flight  is 
examined.  Trim  condition  5  is  a  steady  level  turn  at  a  constant  angular  rate.  Next,  a 
steady  climbing  turn  is  defined  in  trim  condition  6.  Trim  condition  7  is  a  steady 
descending  turn,  the  opposite  of  trim  condition  6.  Trim  condition  8  corresponds  to  a 
barrel  roU.  Next,  trim  condition  9  describes  a  steady  turn  orbit  with  a  weapon  pointing  at 
a  stationary  target.  Finally,  trim  condition  10  defines  a  steady  level  pitch- over. 

Equations  for  each  trim  condition  are  now  developed. 

Trim  Condition  1. 

Trim  condition  I,  steady,  straight  and  level  flight  is  defined  as  y  =  (|)v  =  H  =  0° . 
The  pitch  angle,  0,  is  therefore  equal  to  the  angle  of  attack,  a.  The  (|)v  is  zero  and  the 
heading  angle,  H,  is  equal  to  the  sideslip  angle,  (3 .  Sideshp  angle  is  assumed  to  be  zero 
for  trim  condition  I,  so  H  is  also  zero.  The  Euler  angles  are  therefore  reduced  to 


0  =  a 


4)  =  0° 


(32) 


\\f  =  0° 
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Trim  Condition  2. 


Figure  5.  Aircraft  Climb  Definition 


A  straight  and  steady  climb,  trim  condition  2,  is  obtained  when  (|)v  =  H  =  0°  and  y 
is  specified  greater  than  zero.  From  equations  (26)  and  (27)  0  is  a  function  of  a,  (3  and  y, 
but  (3  is  assumed  to  be  zero.  From  equations  (28),  (29),  (30)  and  (31)  (|)v  is  zero  since  P  is 
zero  and  t|/  is  also  zero.  The  Euler  angles  then  simphfied  to 

sin  0  =  sin  y  cos  a  +  cos  y  sin  a 
cos  0  =  cos  y  cos  a  -  sin  y  sin  a 

4)  =  0°  (33) 

cos  \|/  =  1 
sin  \|/  =  0 
\\t  =  0° 


Trim  Condition  3. 

Trim  condition  3,  a  straight  and  steady  descent,  is  (|)  v  =  H  =  0°  and  y  <  0° .  The 
equations  are  identical  to  trim  condition  2,  (Equations  33),  but  the  quadrant  check  is 
necessary  to  find  the  correct  0  angle. 
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Trim  Condition  4. 


Figure  6.  Aircraft  Crabbing  Flight  Definition 


A  non- zero  sideshp  angle,  crabbing  flight,  is  described  in  trim  condition  4. 
Trimmed  crabbing  flight  is  deflned  by  y  =  H  =  0°  and  (3  0° .  The  angle  (|)  v  is 

determined  from  the  outer  loop  equations  for  the  specified  (3 .  The  Euler  angles  are  found 
from  equations  (26)  -  (31)  and  are  given  by 


sin  0  =  cos  (|)v  sin  a  -  sin  (|)v  cos  a  sin  (3 
cos  0  =  -^Cjj^  -i-C2i^  =  7  (cos  a  cosa)^  -i-  (cos(|)^  cos  a  sin  a  -i-  sin  (|)^  sin  af 
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.  ,  c„  sin  (1)  cos  a 

sin  (|)  =  — ^  - 

cos  e  cos  e 

c„  sin  (I),,  sin  a  sin  a  +  cos  (I),,  cos  a 

cos  ^  -  (34) 

cose  cose 

c,,  cos^cos  a  sin  a  +  sin  (!)„  sin  a 

sin  \|/  =  — —  =  - - 

cos  e  cose 

c,,  cosacosa 

cos  \\f  =  — ^  = - 

cos  e  cose 


Trim  Condition  5. 


Figure  7.  Aircraft  Steady  Turn  Definition 

Trim  condition  5,  a  steady,  level  turn,  is  defined  as  y  =  0°  and  H(t)  =  (ISg  t.  Either 
P  or  (|)v  is  specified  for  the  required  turn,  and  the  unspecified  angle  is  determined  by  the 
outer- loop  equations.  The  general  equations  for  a  turn  are  found  first,  and  then 
simpUfying  assumptions  are  made.  The  Euler  angles  are 


sin  0  =  cos  (|)v  sin  a  -  sin  (|)v  cos  a  sin  P 
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.  ,  c„  sin  (1)  cos  a 

sin  (|)  =  — ^  - 

cos  e  cos  e 


c„  sin  6  sin  a  sin  a  +  cos (!)„  cos  a 

cos  (p  =  — ^  =  — — - — - 

cos  e  cose 


(35) 


sinv|/(t) 

cose 

cos  \|/(t)  = 

cos  e 

where  the  rotation  matrix  elements  are  given  by 

Cjj^  =  cos(05q,t)cosacosa  -costf)^  sin(05,j,t)cosasin  a-sin  (j)^  sin(05q,t)sin  a 
C2i^  =sin(05^t)cosacosa  +  cos(l)^cos(05q,t)cosa  sina  +  sin  (j)^  cos(05>j,t)sin  a 


If  the  aircraft  is  assumed  to  have  no  sideshp  in  the  turn  (P  =0°),  the  equations 
simplify  to  the  following 


sin  0  =  cos  (|)v  sin  a 

cos  0  =  ^cos^'Z+rin^'^rsin^T 


sin(|)  = 


cos  (|)  = 


sin(|), 

cose 


costfi^cosa 


cose 


(35a) 


sin  \|/(t)  = 


cos  \|/(t) 


sin  (03q,t)cos  a  +  sin  (j)^  cos(05>j,t)sin  a 


cose 


cos(03q,t)cos  a  -  sin  (f)^  sin  (05>j,t)sin  a 


cose 
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If  the  aircraft  is  assumed  to  have  no  bank  angle  relative  to  the  velocity  in  the  turn 
((|)v  =  0°),  the  equations  simphfy  to  the  following  for  the  Euler  angles 

0  =  a 

4)  =  0°  (35b) 

\|/(t)  =  05^t  +  a 

Trim  Condition  6. 


Figure  8.  Aircraft  Climbing  Turn  Definition 


Trim  condition  6,  a  steady  climbing  turn,  is  defined  as  y  >  0°  and  H(t)  =  050 1. 

Either  (3  or  (|)v  is  specified  for  the  required  turn,  and  the  unspecified  angle  is  determined 
by  the  outer- loop  equations.  The  general  equations  for  a  turn  are  found  first,  and  then 
simplifying  assumptions  are  made.  The  Euler  angles  are 
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sin  0  =  sin  acos  a  cos  a  -  sin  (f)^  cos  acos  a  sin  a  +  cos(j\,  cos  asin  a 


cos  0  =  -Jcjj^  +  C2i^  (where  H(t)  =  050 1) 


sin  (|)  =  — —  (where  H(t)  =  OSg  t) 
cos  e 


cos  (|)  =  (where  H(t)  =  OSg  t)  (36) 

cos  e 


sin\|/(t)  =  (where  H(t)  =  GSg  t) 

cose 


cos  vi/(t)  = (where  H(t)  =  OSg  t) 
cose 


If  the  aircraft  is  assumed  to  have  no  sideshp  in  the  turn  (|3  =  0°),  the  equations 
simplify  to  the  following 


sin  0  =  sin  acos  a  +  cos(|)^  cos  a  sin  a 


sin  (|) 


sin  (f)^  cos  a 


cose 


(36a) 


sin\|/(t)  = 


cos  a  sin  (05>j,  t  )cos  a  -  [cos  (j)^  sin  asin  (05g  t )  -  sin  (f)^  cos  (05q,  t )]  sin  a 


cose 


If  the  aircraft  is  assumed  to  have  no  bank  angle  relative  to  the  velocity  in  the  turn 
((|)v  =  0°),  the  equations  simplify  to  the  following 
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sin  0  =  sin  acos  a  cos  a  +  cos  asin  a 


.  ,  sin  asin  a 

sin  (p  = - 

cose 


(36b) 


sin\|/(t)  = 

cos  a  sin  (03q,t)cosa  cosa  +  cos(l)^cos(o50t)cosacosa  -coscj)^  sin  asin  (O50t)sin  a 

cose 


Trim  Condition  7. 

Trim  condition  7,  a  steady  descending  turn,  is  defined  as  y  <  0°  and  H(t)  =  050 1. 
The  equations  are  identical  to  trim  condition  6,  (Equations  36),  but  a  quadrant  check  is 
necessary  to  find  the  correct  0,  (p  and  t|/(t)  angles. 

Trim  Condition  8. 
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Trim  condition  8,  a  steady  barrel  roll,  is  defined  as  (3  =  H  =  0°  and  (|)  v(t)  =  05^^  t. 


The  angle,  y,  must  also  be  specified.  The  Euler  angles  are 


sin  0(t)  =  sin  a  cos  a  cos  a  +  cos 


cos  a  sin  a 


cos 


0(t)= 


2  2 

+  Co,  — 


.^(cosacosa)^  -  2  cos  a  cos  a  cos  ^[55^  t)sin  a+  [cos^^  t)sin  asina  +  [sin(o5^  t)sin  aj^ 

•  A.A  C32  sin^t)cosa 

sin  (|)(t)  =  — - 

cos  e  cos  e 


c„  sin  asin  a  +  cos(o5.  tlcos  acos a 
cos  (|)(t)  -  _ 


cose 


cose 


(37) 


c,,  sinm  tlsina 

sin^(t)  = 

cose  cose 


c  cos  a  cos  a  -  cos(05.  t)  sin  a  sin  a 
cos  tl/(t)  -  ^  - 


cose 


cose 
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Trim  Condition  9. 


Figure  10.  AC  - 130  Gunship  Maneuver  Definition 


Trim  condition  9,  a  steady  level  turn  with  a  gun  pointed  out  the  side  of  the  aircraft 
at  a  stationary  target,  is  defined  by  y  =  0°  and  H(t)  =  05,j,  t.  The  aircraft  height  (h)  above 

the  ground  is  a  function  of  V,  n  (load  factor),  and  8  (gun  elevation  angle).  Sideshp  is 
assumed  to  be  zero  ((3  =  0°),  producing  the  same  equations  as  a  level  turn,  (Equations 
35a).  From  Figure  10,  the  radius  (q)  of  the  orbit  is  defined  by  the  aircraft  velocity  and 
load  factor.  The  bank  angle  is  a  function  of  the  load  factor.  The  aircraft  height  above  the 
ground  is  a  function  of  orbit  radius,  bank  angle,  and  weapon  depression  angle. 


V" 


(f)^  =arccos 
h  =  q  tan(4>^  +e) 
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The  three  equations  were  solved  together  to  form  an  expression  for  aircraft  height 
(Equation  38).  Sample  calculations  are  shown  in  Appendix  A.  The  equations  combine  to 
form 


1 

1+  ,  =tane  ^ 
1  - -Jn^  -1  tan e  § 


(38) 


Therefore,  the  required  aircraft  height  above  the  target  is  only  a  function  of 
airspeed,  load  factor,  and  weapon  depression  angle.  If  the  weapon  depression  angle  is 
reduced  to  zero,  the  required  aircraft  height  equation  reduces  to  h  =  V  /g.  To  increase  the 
height  above  the  ground,  the  weapon  depression  angle  must  be  increased.  If  8  is  assumed 
to  be  a  non- zero  positive  number,  the  maximum  height  is  found  by  differentiating 

equation  (38)  and  setting  it  equal  to  zero.  Simphfying  the  equation  results  in  nhmax=  V2  . 
Substituting  the  load  factor  for  maximum  height  produces  equations  (39)  for  the 
maximum  height  above  the  target,  maximum  bank  angle  and  radius. 


h 


max 


1-1-  tane 
1-tane  g 


4* 


'^max 


=  45° 


g 


(39) 
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The  hniax  equation  (39)  is  undefined  at  8  =  45°,  therefore  the  maximum  height 
above  the  target  occurs  as  8  approaches  45°  with  a  bank  angle  of  45°  or  essentially 
pointing  straight  down. 

Trim  Condition  10. 


Figure  1 1 .  Aircraft  Pitch- over  Definition 

Trim  condition  W,  a  level  flight  pitch  over,  is  a  subset  of  trim  condition  ]_ 

2<=  1^  0  °).  This  maneuver  is  envisioned  for  use  at  pointing  a  fixed  tine- of- sight 
weapon  from  the  bottom  of  the  aircraft  at  a  stationary  target.  The  pitch  change  rate,  (De, 
is  a  function  of  aircraft  velocity  and  re.  Therefore  re  =  V/tDe  to  point  at  a  stationary  target 
below  the  aircraft. 

Aircraft  Rotation  Rate  Determination 

With  the  Euler  angles  in  the  navigation  frame,  the  aircraft  angular  change  rates 
resolved  in  the  body  axes  are  then  determined.  Steady  level  flight  2<=  1^  0  °),  is  a 
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subset  of  a  steady  turn  0  °).  The  inner-loop  equations  are  solved  for  a 


steady  turn.  From  the  Euler  equation  (1:  Chapter  2). 


'1 
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the  aircraft  angular  rates  are  obtained 
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'P“ 
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Q 

= 
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R 
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For  the  aircraft  in  a  steady  level  pitch  over,  the  equations  become 
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the  aircraft  angular  rates  are  obtained 
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For  the  aircraft  in  a  steady  roU,  the  equations  become 
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the  aircraft  angular  rates  are  obtained 
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Control  Surface  Deflection  Determination 

The  six- degree- of- freedom  equations  of  motion  are  used  to  find  the  aircraft 
control  surface  deflections  from  the  equations  previously  determined.  Figure  12  shows 


the  definition  of  moments  and  forces. 


Moment  Model  Definition. 


Figure  12.  Three  Dimensional  Aircraft  Moments  and  Forces 

The  thrust  term,  T,  from  the  outer  loop  model  was  divided  into  two  terms,  Ti  and 
T2.  Ti  is  the  thmst  from  the  right  side  engine,  while  Tz  is  the  thmst  from  the  left  side 
engine.  The  three  moment  of  inertia  equations  of  motion  are  (1:  Chapter  2) 


^  (CjR  +  C2P)Q  +  C3 L  +  C4N 
^  C5PR  -  Cg (P"  -  R")  +  C7M 
(CgP  -  C2R)Q  +  C4L  +  CgN 


with  the  following  definitions 


V  yy 

C,  = 


(lyy-IzzKz-I 


IxxL-I 

l^xx  ^zzi^xz  _ 


2  =c. 


^xx^zz  ^xz 


=c  =  c 


pr  "qr 


(43) 

(44) 

(45) 
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XX  zz  xz 

L  =  Lp  +  Lr  +  L.  +  +  Lj 

M  =  M„+M,+M,,+  Mq+M,^-M, 
N  =  Np  +  Np+N,  +  N,^+N,^+N, 

T  =  Tj  +T2 


Roll  Moment  Definition. 

The  roll  moment,  L  ,  is  defined  as  all  moments  about  the  Xp-axis. 
of  the  following  terms 

L  p  =  C  P  ■  q  b  (RoU  moment  due  to  aircraft  roU  rate) 


is  composed 
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L  R  =  C,  R  ■  q  b  (Roll  moment  due  to  aircraft  roU  yaw) 


L  p  =  Cj  a  ■  q  b  (RoU  moment  due  to  aircraft  sidesUp  angle) 

=  C.  a  -q  b  (RoU  moment  due  to  aircraft  aileron  deflection) 

a  a  a 

=  C,  a^  q  b  (RoU  moment  due  to  aircraft  mdder  deflection) 

=  y^sina  ^(T2  -  Tj )  (RoU  moment  due  to  asymmetrical  tmst) 

Pitch  Moment  Definition. 

The  pitch  moment,  M,  is  defined  as  aU  moments  about  the  yn-axis.  M  is 
composed  of  the  foUowing  terms 

Mo  =  q  c  (EquUibrium pitch  moment) 

Ma  =  C^,  a  ■  q  c  (Pitch  moment  due  to  aircraft  angle  of  attack) 

M^<  =  C^^S<q  c  (Pitch  moment  due  to  the  change  in  aircraft  angle  of  attack) 

Mq  =  Q  ■  q  c  (Pitch  moment  due  to  aircraft  pitch  rate) 

a^  -  qc  (Pitch  moment  due  to  aircraft  elevator  deflection) 

M^.  =  (z^-cosa^  +  x^sin  a^^  )(Tj  +  T2 )  (Pitch  moment  due  to  thrust) 

Yaw  Moment  Definition. 

The  yaw  moment,  N,  is  defined  as  all  moments  about  the  Zn-axis.  N  is  composed 
of  the  foUowing  terms 


41 


Np=C„pPqb  (Yaw  moment  due  to  aircraft  roll  rate) 


NR=C„Rqb  (Yaw  moment  due  to  aircraft  roU  yaw) 

Np  =  a  ■  q  b  (Yaw  moment  due  to  aircraft  sideslip  angle) 

N.  =  C„  a  •  q  b  (Yaw  moment  due  to  aircraft  aileron  deflection) 

a  'ag 

Nj  =  C„  a^  ■  q  b  (Yaw  moment  due  to  aircraft  rudder  deflection) 

=  y^cosa^(T2  -  Tj)  (Yaw  moment  due  to  asymmetrical  tmst) 

Control  Surface  Deflection  Equations. 

With  moment  equations  of  motion,  aircraft  trim  conditions  are  determined.  In  the 
most  general  trim  condition,  aU  three  angular  accelerations  (^(^and  ^'^)  were  set  to 
zero.  With  the  appropriate  variable  substitutions,  the  equations  of  motion  (43)- (45)  are 


^0 


Cj  PQ-I-Cj  QR-i-C3qb(C[  P-I-Cj  R-I-Cj  a-i-Cj  a^ -I-Cj  a^)-i-C3y.j.sin  a.j.(T2 -Tj)-I- 
+  C4qb(C„  P  +  C„  R  +  C„,a  +  C„,  a^ +C„  JiJ  +  C4yTCOsaT(T2 -T^) 


(46) 


^0  = 

C„/R  +  (R^  -  P^)  +  C3^(C^^  +  a  +  Q  +  a  J  + 

cosa^.  +Xj  sin  a3.)(Tj  -1-T2) 


(47) 


^^0  = 

C„„PQ  QR+  C9qb(C„^  P+  C„  R  +  C„^a  +  C„^_a,  +  a,)  +  c^y^cos  a^ (T2  -  T; )  + 

-i-C4qb(Cj  P-i-C,  R -i-Cj^a-i-Cj  a^ -I-C,  aj-i- C4y3.  sin a3.(T2 -Tj) 


(48) 
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Solving  equations  (46)  and  (48)  for  5  a  and  5r  results  in  the  following  matrix 


**3 

C3C1,  +‘^4C„, 

-1 

's’ 

^9Cn  +  C4C1 

aa  aa 

.a,. 

“  qb 

.Cr. 

where 

Cp=Cj  PQ+C[  QR +  C3y^sina^(T2 -Tj) +  C4yTCOsa^(T2 -Tj)  + 

+  qb[c3 (C P  +  C,  R  +  Q^a)  +  C4(C„^  P  +  C„  R  +  C„^a)] 

c^=C„  PQ-C„  QR +Coy3.eosa3.(T2 -Tj) +  e4y3.sin  a3.(T2 -Tj)  + 

pq  qr 

+  qb[c9(C„_^P  +  C„  R  +  C„^a)  +  C4 (C,^P  +  R  +  C^^a)] 


which  simplifies  to 


^r(^3^1a  ‘*'^4^113  )  ^p(^9^na  "*"^4^15  ) 

df  d|.  r  df  dr 


(49) 


qb[(c3Ci  +C4C  )(C9C  +C4C,  )-(c3Ci  +C4C  )(C9C  +C4C1  )] 


Solving  equation  (47)  for  5e 


6e= 

C7(z3,cosa3. +X3.sina3.)(Tj +T2)-C^  PR-C^^  (R^-P^) 

c^qc 

dg 


(50) 


If  it  is  assumed  that  the  thmst  is  the  same  from  both  engines  (Ti=T2),  Equation 
(50)  simphfies  to  the  lateral- direction  equations 
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(51) 


(^9^na_  ^  ^4^1,  ) 


(^9^nii  “*“^4^13  )  (^3^1a  "*'^4Cna  ) 


qb^CjC,  +C4C  )(C9C  +C4C,  )-(c3C,  +C4C  )(C9C  +C4C,  )] 


where 


Cp  =C,  PQ  +  Q  QR  +  qb[c3(CiP  +  C,  R  +  C,a)+C4(C„P  +  C„R  +  C„a)] 

roT  pq  qr  pra  pra 

c,  =C„  PQ-C„  QR  +  qb[c9(C„P+C„R  +  C„,a)+C4(C,P  +  QR  +  C,,a)] 


The  longitudinal  equation  (50)  is  not  simplified  by  this  assumption,  since  both  engines 
thmst  produce  the  same  longitudinal  moment. 

Summary 

In  this  chapter,  inner- loop  equations  were  developed  to  solve  the  control  surface 
deflections,  from  a  given  set  of  control  variables  and  state  variables.  For  a  desired 
trajectory  and  orientation  the  state  variables  were  determined  and  the  outer- loop 
equations  were  used  to  solve  the  required  control  variables.  The  equations  developed  in 
this  chapter  used  the  values  from  the  outer- loop  and  the  desired  angular  acceleration 
values  in  an  inner- loop  system  to  determine  the  required  control  surface  deflections. 

Therefore,  by  choosing  the  variables  h^ ,  a ,  ,  ((j)^  or  a  ),  H,  05^ ,  GSg  ,  05^  and  Ti  to 

describe  the  flight  path,  the  required  aircraft  settings  are  found.  Now  that  the  equations 
were  developed,  specific  flight  paths  and  orientations  are  examined  and  control  surface 
deflections  calculated  in  the  next  chapter. 
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rv.  Aircraft  Flight  Trajectories 


Overview 

In  Chapter  2,  the  point  mass  model  was  developed  to  provide  a  basis  for  the 
control  laws  for  the  trajectory  of  the  aircraft.  The  outer- loop  equations  were  found  from 
the  nonlinear  equations  of  motion  and  apphed  to  the  actual  Learjet-25  parameters.  In 
Chapter  3,  the  outer- loop  equation  outputs  were  used  with  the  transformation  matrices 
and  the  remainder  of  the  equations  of  motion  to  find  the  control  surface  deflections 
required  to  achieve  the  desired  fiight  trajectory  and  orientation.  In  Chapter  4,  the 
equations  developed  in  the  previous  two  chapters  are  used  to  simulate  selected  fiight 
trajectories.  By  choosing  the  correct  input  states,  the  required  control  surface  deflections 
for  the  desired  fiight  path  are  found.  Selected  aircraft  fiight  trajectories  are  examined 
using  the  control  laws  developed  in  the  last  two  chapters.  By  choosing  the  correct  input 
states,  the  desired  fiight  path  and  orientation  for  the  Learjet-25  is  determined.  The 
impacts  of  variations  in  the  stabihty  derivatives  are  investigated  in  Chapter  5  to 
determine  which  parameters  require  the  most  accuracy.  Finally,  in  Chapter  6  the 
theoretical  values  tested  using  the  Learjet-25  in  flight  tests  are  compared  to  the 
theoretical  calculations. 

The  control  laws  are  developed  to  provide  an  automatic  control  system  for  the 
aircraft.  The  goal  was  to  provide  simple  inputs,  to  produce  the  necessary  performance 
from  the  aircraft.  A  Matlab®  routine  was  developed,  using  the  aircraft  control  laws,  to 
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solve  the  equations  from  the  state  variable  inputs.  The  Matlab®  commands  are  shown  in 
the  Appendix  B. 

Trim  Condition  Investigation 

In  this  chapter,  input  state  variable  scenarios  are  examined  to  determine  the 
predicted  aircraft  behavior.  The  results  are  compared  to  the  equations  from  Chapter  3. 
Each  of  the  trim  conditions  is  solved  for  the  required  control  surface  deflections.  The 
most  basic  case  examined  is  trimmed  steady,  level  flight  on  a  straight  trajectory,  trim 
condition  ]_.  Next,  trim  conditions  2  and  3,  climbs  and  descents,  are  explored  by  varying 
the  aircraft  flight  path  angle.  Aircraft  crabbing  (trim  condition  4),  flight  with  non- zero 
sideshp,  is  investigated  by  varying  |3  and  the  resulting  (|)v.  Then  the  different  definitions 
of  a  steady- level  turn  are  examined  by  solving  the  variations  of  trim  condition  5.  Next, 
trim  conditions  6  and  7,  trimmed  climbing  and  descending  turns,  are  examined.  The 
barrel  roU,  trim  condition  8,  is  examined  by  setting  the  change  in  roll  angle  equal  to  a 
constant.  Next,  the  special  case  of  the  AC  - 130  gunship  flight  maneuver  (trim  condition 
9)  is  studied.  The  flight  trajectory  is  defined  by  a  steady,  level  turn  with  the  gun  barrel 
pointing  out  the  side  of  the  aircraft  pointing  at  a  stationary  target.  Trim  condition  10,  a 
constant  rate  pitch  over,  is  the  last  trimmed  flight  maneuver  evaluated.  Each  trim 
condition  maneuver  are  now  detailed. 

Trim  Condition  1  Investigation. 

The  first  flight  condition  examined  was  the  most  basic  trimmed  condition:  level 
flight  with  no  turning.  AH  the  input  states,  a,  (j)^,  a,  (^5  ^fand  are  set  equal  to  zero. 
Eirst  the  vahdity  of  the  equations  over  the  operational  envelope  was  investigated.  After 
finding  the  trimmed  steady,  level  flight  condition,  the  remaining  flight  input  parameters 
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were  varied  to  show  the  affect  on  the  control  surface  deflection  values.  The  altitude  was 


varied  between  9,000  and  20,000  feet  (ft),  to  determine  the  effect  on  the  trimmed  flight 
conditioa  By  changing  the  altitude,  the  standard  atmosphere  density  and  temperature 
changes,  which  affects  the  amount  of  elevator  deflection  (5e)  and  thrust  (T)  required  from 
the  engines.  Next,  the  influence  of  aircraft  weight  on  T  and  5e  was  investigated  by 
varying  the  aircraft  gross  weight.  The  airspeed  effects  were  then  examined  by  varying 
the  aircraft  velocity  from  the  equihbrium  value.  Finally,  asymmetrical  thmst  was 
analyzed  by  varying  the  amount  of  thmst  from  the  right  engine  at  the  trimmed  steady 
level- flight  condition.  Asymmetrical  thmst  could  be  used  as  a  control  input  to  create 
yawing  moments  about  the  center  of  gravity. 

A  center  test  point  flight  condition  was  chosen  to  standardize  the  flight 
parameters.  A  standard  aerospace  atmosphere  and  constant  gravitational  acceleration 
were  chosen  to  perform  all  the  calculations  since  simplified  models  exist  and  are 
typically  used  to  allow  comparisons  between  test  conditions.  The  aircraft  derivatives 
were  assumed  valid  from  9,000  to  20,000  ft,  the  standard  test  altitude  was  chosen  as 
15,000  ft.  Standard  atmosphere  equations  were  used  to  calculate  the  corresponding 
density  and  temperature  for  the  selected  altitude.  Also,  the  derivatives  were  assumed 
apphcable  for  aircraft  velocities  between  237  knots  indicated  airspeed  (KIAS)  and  326 
KIAS.  Therefore  267  KIAS  (approximately  450  feet  per  second)  was  chosen  as  the 
standard  test  velocity.  The  center  test  point  operating  weight  for  the  Leaijet-25  was 
chosen  to  be  13,500  pounds  (lbs). 

The  control  law  equations  were  solved  at  the  center  test  point  conditions.  The 
mdder,  5r,  and  aileron  deflections,  5  a,  are  zero,  since  the  aircraft  does  not  produce  lateral- 
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directional  moments  in  steady  level  flight.  The  required  elevator  deflection,  6e,  was  -0.02 
radians  (rad)  or  - 1 . 17  degrees  (deg),  for  level  flight  at  the  center  test  point  condition.  It 
was  found  that  the  aircraft  nose  must  point  up  shghtly  to  maintain  level  flight,  so  the 
angle  of  attack,  a,  and  the  pitch  angle,  0,  were  both  0.0782  rad  (4.48°).  The  angles  a  and 
0  are  the  same  since  all  other  aircraft  angles  are  zero.  These  values  match  the  equations 
(32)  found  in  Chapter  3.  For  steady  level  flight  at  the  standard  test  condition,  the  aircraft 
required  1,053  lbs  of  thrust.  The  angular  velocities  Q,  and  R )  and 

accelerations  (^^and  it.^),  and  aircraft  angles  (a,  a,(l)^ ,  and  H  )  were  all  zero.  This 
flight  condition  is  used  as  the  baseline  for  comparison  for  the  remainder  of  the  flight 
scenarios  examined.  The  first  flight  parameter  explored  is  the  altitude  effect. 

Altitude  Effect  Investigation. 


verses  Altitude  (h) 
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The  control  law  equations  were  solved  for  various  altitudes  over  the  apphcable 
envelope.  Changing  the  altitude  affects  the  required  thrust  of  the  aircraft.  The  higher  the 
altitude,  the  lower  the  temperature  and  air  density,  which  lowers  the  required  thmst. 
Therefore,  the  only  parameters  that  changed  were  a,  0,  T,  and  the  required  5e.  0  was 
equal  to  a  at  all  points.  Figure  13  shows  the  relationship  between  altitude  and  the 
required  5e  and  the  aircraft  angle  of  attack,  a.  The  required  elevator  deflection  varied  up 
to  2.38°  over  the  range  of  vahd  altitudes.  A  nearly  hnear  relationship  was  observed,  with 
-0.22°  of  5e  required  for  every  1,000  ft  increase  in  altitude.  A  similar  relationship  was 
observed  for  a,  with  a  0.21°  increase  in  a,  for  every  1,000  ft  increase  in  altitude.  The 
required  thrust  was  also  varied  over  the  altitude  envelope. 
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The  linear  relationship  between  altitude  and  required  thrust  is  shown  in  Figurel4. 
Changing  altitude  had  a  moderate  effect  on  the  thmst,  with  a  350  lbs  difference  between 
the  top  and  the  bottom  of  the  altitude  envelope.  For  each  1,000  ft  increase  in  altitude, 
32.4  lbs  less  of  thmst  was  required  to  maintain  trimmed  flight.  This  result  was  expected 
since  the  lower  density  at  a  higher  altitude  reduces  the  drag,  which  reduces  the  required 
thmst. 

Weight  Effect  Investigation. 


Figure  15.  Elevator  Control  Surface  Deflection  (5e)  and  Angle  of  Attack  (a) 

verses  Aircraft  Weight  (W) 


After  examining  the  effect  of  altitude  on  the  trim  condition,  the  aircraft  gross 
weight  was  varied  around  the  center  weight  of  13,500  lbs  to  find  the  impact  on  the 
control  settings.  The  aircraft  weight  range  was  chosen  between  12,000  and  15,000  lbs. 
The  results  are  shown  in  Figure  15.  Changing  the  aircraft  weight  changed  the  amount  of 
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lift  required,  which  was  directly  related  to  a.  The  elevator  angle  was  varied  by  changing 
a.  Therefore,  the  heavier  the  aircraft,  the  more  negative  5e  and  positive  a  required  for 
trimmed  level  flight.  The  slope  of  the  weight  tine  was  greater  than  for  the  altitude  study: 
a  1,000  lbs  increase  in  weight  requires  a  0.34°  decrease  in  the  elevator  deflection  angle. 
A  similar  relationship  was  found  in  the  a,  a  1,000  lbs  increase  in  weight  requires  a  0.32° 
increase  in  angle  of  attack. 


Figure  16.  Thrust  (T)  verses  Aircraft  Weight  (W) 

Again  a  hnear  relationship  was  noted  for  thmst  over  the  weight  flight  envelope 
(see  Figure  16).  Weight  changes  have  a  limited  impact  on  the  required  thmst  of  the 
Learjet-25.  A  1,000  lbs  increase  in  weight  results  in  only  a  23.1  lbs  increase  in  thmst. 
The  variation  was  caused  by  the  increased  weight,  which  increased  the  required  a  to 
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produce  more  lift,  which  increased  drag.  Now  that  the  effeet  of  ehanging  the  aircraft 
weight  is  known,  ehanges  in  veloeity  are  presented. 

Velocity  Effect  Investigation. 


Figure  17.  Elevator  Control  Surface  Deflection  (5e)  and  Angle  of  Attack  (a) 

verses  Aircraft  Velocity  (V) 


The  next  steady  level  flight  parameter  examined  was  the  affect  of  a  change  in 
velocity  from  the  equilibrium  veloeity  of  450  feet  per  seeond  (ft/s).  The  test  points  were 
chosen  over  the  assumed  vahd  range  for  the  aircraft  derivatives,  from  237  KIAS  (400 
ft/s)  to  326  KIAS  (550  ft/s).  To  aehieve  the  velocity  change,  while  maintaining  steady 
flight,  the  trimmed  thmst  setting  was  ehanged  to  vary  velocity.  This  resulted  in  a  change 
in  a,  and  the  required  5e.  The  rest  of  the  eontrol  variables  were  not  affected.  The  results 
are  shown  in  Figure  17.  Second-order  relationships  with  velocity  were  seen  for  both  5e 
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and  a.  As  the  velocity  was  changed  from  400  ft/s  to  the  upper  limit  of  550  ft/s,  the 
required  5e  was  increased  2.89°  and  a  was  decreased  2.73°.  Therefore,  increasing  the 
airspeed  required  a  smaller  negative  5e  to  maintain  trimmed  level  flight,  with  a 
corresponding  decrease  in  a. 


Changing  the  velocity  from  the  steady  state  value  had  a  moderate  impact  on  the 
required  aircraft  thmst,  as  shown  in  Figure  18.  The  parameters  have  a  linear  relationship 
with  a  slope  of  3.2  lbs  of  thmst  per  ft/s.  When  increasing  the  airspeed  from  the  lower 
limit  of  400  ft/s  to  the  upper  limit  of  550  ft/s,  478  lbs  more  thmst  was  required.  The 
result  was  expected,  since  higher  velocity  results  in  more  drag  requiring  more  thmst.  The 
first  three  parameters  affected  the  longitudinal  stabflity  of  the  aircraft;  next  a  lateral- 
direction  parameter  will  be  presented. 
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Asymmetrical  Thrust  Effect  Investigation. 


Figure  19.  Aileron  Surface  Deflection  (5a)  and  Rudder  Surface  Deflection  (5r) 
verses  Percent  of  Right  Engine  Thmst 

The  Learjet  has  two  engines,  on  each  side  of  the  center  of  gravity.  If  the  thmst 
from  each  engine  is  not  the  same,  roll  and  yaw  moments  wiU  be  created  about  the  center 
of  gravity.  Asymmetrical  thmst  was  examined  by  calculating  the  percentage  of  thmst  out 
of  the  right  engine.  An  input  of  100%,  represents  aU  the  thmst  was  from  the  right  engine, 
with  0%  from  the  left  engine.  For  most  calculations  in  this  study  it  was  assumed  that 
50%  of  the  tmst  was  from  each  engine.  The  pitch  control  variables  (5e  and  a)  were  not 
changed  with  a  thmst  asymmetry  since  the  longitudinal  forces  and  moments  remained 
constant.  To  maintain  steady  level  flight,  with  asymmetrical  thmst,  the  aileron  (5  a)  and 
mdder  (5r)  were  deflected.  Figure  19  reproduces  the  relationship  between  asymmetrical 
thmst  and  the  required  lateral  control  surface  deflections.  Finear  relationships  were 
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observed  for  both  control  surfaces.  For  a  one  percent  increase  in  thrust  from  the  right 
engine  over  the  left  engine,  a  0.01°  decrease  in  5a  and  a  0.05°  decrease  in  5r  were 
required.  Thus,  if  60%  of  the  total  thmst  was  from  the  right  engine  (a  10%  increase),  5a 
must  be  decreased  0.13°  and  5r  decreased  0.53°  from  the  symmetrical  test  conditioa 
Now  that  steady  level  flight  has  been  examined,  some  other  simple  trimmed  flight 
maneuvers  will  be  studied. 

Trim  Conditions  2  and  3  Investigation. 

The  next  aircraft  maneuvers  examined  were  straight  and  steady  trimmed  climbs 
and  descents.  The  aircraft  was  assumed  to  remain  in  the  same  orientation,  so  the  angular 
rotations  and  accelerations,  P,  Q,  R,  ^fand  were  all  still  zero.  Initially  the 

maneuver  was  limited  to  the  longitudinal  plane,  so  that  a,  (|)^ ,  and  H  were  also  zero.  A 

steady  change  in  altitude  was  achieved  by  setting  the  flight  path  angle,  y,  equal  to  a  non¬ 
zero  value  and  solving  the  control  law  equations  for  trimmed  flight. 


Figure  20.  Elevator  Surface  Deflection  (5e),  Angle  of  Attack  (a)  and  Pitch  Angle 


(0)  verses  Flight  Path  Angle  (y) 


55 


From  steady  straight  and  level  flight,  the  aircraft  maneuver  examined  was  a 
change  in  altitude  by  a  trimmed  climb  or  descent.  The  input  parameter,  y,  was  varied 
while  the  rest  of  the  input  parameters  were  unchanged.  The  climb  angle  was  varied  until 
the  upper  and  lower  limits  of  the  available  thrust  were  found.  The  required  angle  of 
attack  and  pitch  angle  were  also  calculated.  The  maneuver  was  limited  to  the 
longitudinal  plane;  therefore  5a  and  5r  remained  zero.  The  graph  of  0,  5e  and  a  are  shown 
in  Figure  20.  Angle  of  attack  and  elevator  deflection  remain  nearly  constant  for  aU  y. 

The  pitch  angle  was  linearly  related  to  y  with  a  1.0°  change  in  y  resulting  in  the  same 
change  in  0. 


Figure  21.  Thrust  (T)  verses  Flight  Path  Angle  (y) 

The  rate  of  climb  is  directly  related  to  the  amount  of  required  thmst.  Flight  path 
angle  and  thmst  are  shown  in  Figure  21.  Gamma  and  thmst  are  a  non-linear  function  of 
sine  and  cosine,  which  flatten  out  near  y=0°,  and  increase  as  y  increases.  Small  changes 
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in  flight  path  angle  produce  a  large  change  in  required  thrust,  especially  the  further  from 
level  flight.  The  rate  of  climb,  and  associated  flight  path  angle,  require  more  thmst  as  y 
increases.  Also  for  a  trimmed  descent,  httle  to  no  thmst  is  required.  The  lower  limit  of  y 
corresponds  to  zero  thmst  required  to  remain  trimmed.  Next  a  lateral- directional  trim 
condition  was  investigated,  during  crabbing  flight. 

Trim  Condition  4  Investigation. 

Level  crabbing  flight  with  a  non- zero  sideshp  angle  was  examined  in  trim 
condition  4.  The  aircraft  was  assumed  to  be  in  level  flight,  with  a  constant  heading.  The 

forces  and  moments  were  balanced,  therefore  P,  Q,  R,  §5and  l¥  were  all  zero. 

The  angles  y  and  H  were  also  set  equal  to  zero.  The  value  of  (3  began  at  zero  and  was 
increased  in  magnitude  in  both  the  positive  and  negative  directions  until  a  limit  was 
reached.  Since  (3  was  specified,  (|)v  was  determined  from  the  outer- loop  equations  to 
balance  the  forces  and  moments.  The  Learjet-25  had  a  ±10°  (3  limit,  which  was  the 
limiting  factor  in  the  sideshp  inveshgation.  The  longitudinal  parameters  T,  a,  0,  and  5e 
were  essentially  constant  for  the  entire  range  of  |3  examined. 
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Sideslip  Angle  Effects 


+&■ 


—I - 1-10  I  '  ^ 

Sideslip  Angie  (deg.) 


Figure  22.  Aileron  Surface  Deflection  (5a),  Rudder  Surface  Deflection  (5r)  and 
Bank  Angle  ((|))  verses  Sideslip  Angle  ((3) 


In  straight  and  level  flight,  (3  was  varied  to  the  limits,  while  the  rest  of  the  input 
variables  were  held  constant.  The  lateral- directional  parameters  were  calculated  and 
plotted  in  Figure  22.  Aileron  deflection  was  inversely  proportional  to  the  sideshp  angle, 
with  a  0.79°  decrease  in  5a  for  every  1.0°  increase  in  |3.  Conversely,  5r  increased  as  (3 
increased,  with  a  1.63°  increase  in  5r  for  each  1.0°  increase  in  sideshp  angle.  FinaUy, 
bank  angle  was  linearly  related  to  [3 .  A  0.70°  increase  in  (|)  corresponded  to  a  1 .0° 
increase  in  (3 .  Next,  another  lateral- directional  trim  condition,  a  steady,  level  turn,  was 
investigated. 

Trim  Condition  5  Investigation. 

Level  turns  were  investigated  by  varying  bank  angle  and  sideshp  angle.  The 
aircraft  was  assumed  to  be  in  level  flight  with  a  constantly  changing  heading,  H  =  050 1,  a 
function  of  time.  The  (|)v  and  (3  were  determined  by  specifying  one  variable  and 
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calculating  the  other  by  the  outer- loop  equations.  The  control  surface  deflections  are 
constant  with  respect  to  time.  In  Figure  23,  (|)v  was  specified.  For  trim  condition  5,  y  was 
also  set  equal  to  zero.  The  parameters  T,  a,  0  and  6e  did  not  change  significantly  when  (|)v 
was  varied.  The  limiting  factor  was  the  maximum  5r.  The  bank  angle  was  varied  both 
directions  until  the  mdder  limit  was  exceeded. 


Figure  23.  Level  Flight  Aileron  Surface  Deflection  (5  a),  Rudder  Surface 
Deflection  (5r)  and  Sideshp  Angle  ((3)  vs.  Bank  Angle  ((|)v) 


From  Figure  23,  there  was  a  linear  relationship  between  (|)  v  and  5a,  5r,  and  [3 . 
Aileron  deflection  was  inversely  proportional  to  (|)v,  with  a  1.0°  increase  in  bank  angle 
decreasing  5a  by  1.38°.  On  the  other  hand,  when  (|)v  was  increased  1.0°,  5r  increased 
2.85°.  Finally,  the  slope  of  the  [3  per  (|)v  line  was  1.75°/°.  Next  the  more  general  cases  of 
steady  turns  were  investigated,  when  the  flight  path  angle  was  not  zero. 
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Trim  Condition  6  Investigation. 

For  trim  condition  6,  a  turn  rate  and  climb  angle  were  chosen  to  explore  the  effect 
of  bank  angle  on  a  climbing  turn.  As  in  trim  condition  5,  the  heading  angle  was 
described  by  H  =  Cdg  t.  Additionally,  y  was  defined  as  a  positive  angle.  As  before,  (|)v 

was  varied  until  an  aircraft  parameter  limit  was  reached,  and  [3  was  calculated  from  the 
outer- loop  equations.  The  mdder  deflection  was  still  the  limiting  factor  for  this  case. 

The  pitch  angle  changed  shghtly  over  the  (|)v  range  examined,  so  0  was  also  plotted. 


Figure  24.  Climbing  Flight  Aileron  Surface  Deflection  (5  a),  Rudder  Surface 


Deflection  (5r)  and  Sideshp  Angle  ((3)  vs.  Bank  Angle  ((|)v) 


The  aircraft  parameters  for  trim  condition  6  were  nearly  identical  to  trim 
condition  5.  The  data  were  plotted  in  Figure  24.  The  slope  of  5a/(|)v  and  (3/(|)v  did  not 
change.  The  slope  of  5i/(|)v  decreased  from  2.85°/°  to  2.84°/°  for  this  case.  The  pitch 
angle  decreased  shghtly  as  the  magnitude  of  (|)v  increased.  The  0  began  at  7.84°  and 
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decayed  to  3.58°  at  (|)v  =±1 1 .46° .  Next  a  descending  turn  was  investigated  by  restricting  y 


to  negative  values.  Thrust  results  were  nearly  identical  to  the  flight  path  angle 


investigation,  as  expected. 


Trim  Condition  7  Investigation. 


Figure  25.  Descending  Flight  Aileron  Surface  Deflection  (5a),  Rudder  Surface 


Deflection  (5r)  and  Sideslip  Angle  ((3)  vs.  Bank  Angle  ((|)v) 


The  next  aircraft  maneuver  investigated  was  trimmed  descending  turns.  The 
assumptions  and  limitations  for  trim  condition  7  are  identical  to  trim  condition  6  with  one 
exception:  y  was  negative.  The  results  are  shown  in  Figure  25.  The  slopes  of  all  three 
linear  relationships  were  unchanged  from  trim  condition  6.  The  primary  difference 
between  this  case  and  the  previous  trim  condition  was  the  effect  of  y  on  0.  Lowering  the 
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y,  shifted  the  0  curve  down  and  flattened  the  curve,  making  0  less  of  a  function  of  y.  The 


next  trim  condition  examined  was  the  rolling  effects  during  a  barrel  roU. 

Trim  Condition  8  Investigation. 


Barrel  Roll  Effects 


Figure  26.  Barrel  Roll  Aileron  Surface  Deflection  (5a)  and  Rudder  Surfaee 
Deflection  (5r)  vs.  Bank  Angle  ((|)v) 


A  barrel  roll,  with  eonstant  roll  rate  about  an  axis  above  the  Xb-axis,  was  modeled 
in  trim  eondition  8.  The  barrel  roll  was  defined  as  (3  =  H  =  0°  and  (|)v(t)  =  05^  t.  The 

flight  path  angle  was  also  set  to  a  positive  number.  By  restrieting  (3  to  zero,  the  radius 
from  the  axis  of  rotation,  rb,  was  zero.  Therefore  trim  condition  8  was  the  roU  rate  of  the 
aircraft.  To  increase  rb,  (3  must  be  increased  to  a  positive  number.  The  larger  (3 ,  the 
larger  rb.  Non- zero  |3  was  not  investigated.  The  parameters  a,  0,  T,  and  5e  remained 
constant.  The  lateral- directional  control  surface  deflections  were  linear  functions  of  03^^ 
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(see  Figure  26).  The  roll  rate-limiting  factor  was  the  amount  of  5a  available.  The 
influence  of  roU  rate  on  5r  was  small;  a  0.01  deg/s  increase  in  05^  created  a  0.64° 

decrease  in  5r.  Meanwhile,  a  0.01  deg/s  increase  in  05^^  produced  a  3.39°  decrease  in  5a. 

The  Learjet-25  cannot  complete  a  roU.  This  investigation  only  examined  the  trimmed 
condition  to  start  the  roU  from  wings  level.  While  a  complete  roU  may  be  mathematicaUy 
possible,  the  maneuver  was  not  physicaUy  possible  with  this  aircraft  model.  In  the  next 
trim  condition,  an  operationaUy  representative  AC  - 130  gunship  orbit  was  studied. 

Trim  Condition  9  Investigation. 

The  first  mission  representative  aircraft  maneuver  examined  was  the  orbit  of  the 
AC- 130  gunship.  The  orbit  trajectory  was  described  by  a  trimmed  level  steady  turn  (trim 
condition  5).  The  parameters  5a,  5r,  and  (3  are  aU  functions  of  (|)v.  Figure  23  also  appUes 
to  trim  condition  9.  For  the  AC  - 130  orbit,  the  aircraft  height  above  the  target,  orbit 
radius  and  bank  angle  were  varied  as  functions  of  load  factor  and  weapon  depression 
angle.  The  velocity  was  held  constant  at  450  ft/s.  Figure  27  shows  the  effect  on  height 
above  the  target,  and  Figure  28  shows  the  relationship  of  p  and  (|)v  to  n. 
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Figure  27.  AC- 130  Gunship  Weapon  Depression  Angle  (£)  and  Load  Factor  (n) 

verses  Height  (h) 

The  height  of  the  gunship  orbit  above  the  target  is  a  function  of  weapon 
depression  angle  and  load  factor  (see  Figure  27).  For  the  weapon  depression  angle  curve, 
the  load  factor  was  fixed  at  1.15,  which  resulted  in  an  orbit  radius  of  11,083  ft  and  a  bank 
angle  of  30°.  As  the  weapon  depression  angle  asymptotically  approached  45°  the  height 
above  the  target  approached  infinity.  For  the  specified  fiight  condition,  the  minimum 
height  above  the  target  was  6,294  ft.  For  the  load  factor  curve,  the  weapon  depression 
angle  was  set  at  20°,  and  n  was  varied  from  1.0  to  1.4  (the  maximum  bank  angle).  The 
height  approached  a  minimum  at  12,850  ft  and  increased  slightly  as  n  approached  1.4. 

As  the  load  factor  approached  1.0,  the  height  above  the  target  approached  infinity. 

Utilizing  the  load  factor  curve,  the  weapon  depression  angle  curve  would  move  shghtiy 
right  as  load  factor  increased  and  dramatically  right  as  n  approached  1.0  since  n=1.15  was 
near  the  minimum  of  the  load  factor  curve. 
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Figure  28.  AC- 130  Gunship  Orbit  Bank  Angle  ((|)v)  and  Radius  (p)  verses 

Height  (h) 


The  effeet  of  load  faetor  on  orbit  radius  and  bank  angle  is  shown  in  Figure  28. 
The  orbit  radius  curve  asymptotically  approached  5,000  ft  as  load  factor  approached  the 
upper  limit.  Also,  as  n  approached  1.0,  the  orbit  radius  approached  infinity  (as  expected 
since  n=1.0  is  level  fiight).  The  bank  angle  curve  approached  0.0°  as  n  approached  1.0. 
Likewise,  at  the  upper  limit  as  n  approached  1.4,  (|)v  approached  45°.  Both  curves  are 
independent  of  8  and  apply  to  aU  trim  condition  9  flight  conditions. 
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Trim  Condition  10  Investigation. 


Figure  29.  Pitch- Over  Elevator  Surface  Deflection  (5e)  verses  Pitch  Rate  (Q) 

The  flnal  aircraft  maneuver  investigated  was  the  trimmed  pitch- over.  The 
purpose  of  trim  condition  1^  was  to  create  a  maneuver  where  a  device  mounted  to  the 
bottom  of  the  aircraft  was  pointed  at  a  fixed  target  below  the  aircraft.  By  deflning  the 
correct  V  and  03. ,  a  radius,  re,  was  defined  for  the  trajectory  arc.  The  pitch  rate  dictates 
the  amount  of  time  the  device  points  at  the  target.  A  slower  pitch  rate  produced  a  longer 
time  on  target.  The  maneuver  was  defined  as  level  with  '^0°.  Level  flight  with  Y=0°, 
produced  a  constantly  changing  a.  For  a  slow  pitch  rate,  level  flight  would  not  keep  the 
device  pointed  at  a  stationary  target.  The  pitch  angle  varied  continuously  throughout  the 
maneuver.  The  maneuver  was  also  limited  to  the  longitudinal  plane;  therefore  all  lateral- 
directional  parameters  were  neghgible.  The  primary  parameter  of  interest  for  trim 
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condition  W  was  5e.  The  data  were  plotted  in  Figure  29.  The  pitch  rate  was  varied  from 
0.0  deg/s  to  the  elevator  deflection  limit  15  deg/s.  The  elevator  control  surface  deflection 
was  a  hnear  function  of  the  pitch  rate.  A  0.01  deg/s  increase  in  pitch  rate  resulted  in  a 
10.27°  decrease  in  the  required  5e. 

Summary 

Individual  trim  conditions  were  investigated  in  Chapter  4.  For  straight  and  level 
flight,  nearly  linear  relationship  were  found  for  5e  and  thmst  for  both  altitude  and  weight 
changes.  For  velocity  changes,  a  hnear  relationship  was  found  between  velocity  and  the 
thmst  setting.  Velocity  changes  produced  a  second  order  relationship  with  5e  though. 

The  asymmetric  thmst  investigation  found  linear  functions  for  both  5a  and  5r.  Climbs  and 
descents  were  found  to  have  a  direct  impact  on  thmst  setting,  while  a  negligible  impact 
on  5e.  5a  and  5r  were  linearly  changed  by  changes  in  the  sideslip  angle.  For  steady  turns, 
a  linear  relationship  was  observed  between  5a  and  5r  with  changes  in  bank  angle.  The 
same  relationship  was  observed  between  5a  and  5r  with  bank  angle  as  the  flight  path  angle 
was  varied.  For  a  barrel  roU,  linear  relationships  were  found  for  both  5  a  and  5r  as  roll  rate 
was  varied.  The  effects  of  load  factor  and  weapon  depression  angle  on  height  above  the 
target  were  determined  for  a  simulated  AC  - 130  gunship  orbit.  Also  the  impacts  of  bank 
angle  and  orbit  radius  on  load  factor  were  plotted.  Finally  a  hnear  relationship  was 
observed  between  5e  and  pitch  rate  for  the  pitch-over  maneuver. 
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V.  Stability  Derivative  Sensitivity  Analysis 


Overview 

In  Chapter  2,  the  point  mass  model  was  developed  to  provide  a  basis  for  the 
control  laws  for  the  trajectory  of  the  aircraft.  The  outer- loop  equations  were  found  from 
the  nonlinear  equations  of  motion  and  apphed  to  the  actual  Learjet-25  parameters.  In 
Chapter  3,  the  outer- loop  equation  outputs  were  used  with  the  transformation  matrices 
and  the  remainder  of  the  equations  of  motion  to  find  the  control  surface  deflections 
required  to  achieve  the  desired  fiight  trajectory  and  orientation.  In  Chapter  4,  the 
equations  developed  in  the  previous  two  chapters  were  used  to  simulate  selected  fiight 
trajectories.  By  choosing  the  correct  input  states,  the  required  control  surface  deflections 
for  the  desired  flight  path  were  found.  The  impacts  of  variations  in  the  stabihty 
derivatives  are  investigated  in  Chapter  5  to  determine  which  parameters  require  the  most 
accuracy.  Each  stabihty  derivative  is  varied  ±40%  to  determine  the  change  in  predicted 
control  surface  deflection  value.  Control  surface  deflection  changes  greater  than  1.0° 
from  the  original  value  are  deemed  significant.  The  sensitivity  analysis  was  used  to 
identify  the  most  critical  parameters  and  possible  sources  of  error.  FinaUy,  in  Chapter  6 
the  theoretical  values  tested  using  the  Leaijet-25  in  flight  tests  are  compared  to  the 
theoretical  calculations. 

Sensitivity  Analysis 

The  stabihty  derivatives  were  divided  into  two  groups:  longitudinal  stabihty 
derivatives  and  lateral- directional  stabihty  derivatives.  The  longitudinal  stabihty 
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derivatives  were:  ,  and 

The  longitudinal  stability  derivatives  were  investigated  in  two  trimmed  flight  conditions 
(straight  and  level  flight  and  climbing  flight).  The  lateral- directional  stabihty  derivatives 


were 


•  C  C  C  C  C  0000  0  0  0  00  and 

•  Ya  ’  yp  ’  Yr  ’  ^Yaa  ’  ^  Yar  ’  ’  ‘p  ’  4  ’  ^  W  ’  ^  lar  ’  ^  n  > 


.  The  lateral- directional  stability  derivatives  were  investigated  in  two  trimmed  flight 
conditions  (crabbing  flight  and  steady- turn  flight).  Straight  and  level  flight  did  not 
produce  any  forces  in  the  lateral- directional  axes,  and  no  change  in  the  control  surface 
deflections. 

Longitudinal  Stability  Derivative  Sensitivity. 

The  longitudinal  stabihty  derivatives  sensitivities  were  found  for  two  cases.  Each 
stabihty  derivative  was  varied  ±40%  from  the  original  value  to  estimate  the  effect  of 
errors  in  the  stabihty  derivatives  on  the  predicted  control  surface  deflections.  Stabihty 
derivatives  calculated  from  wind  tunnel  data  may  have  up  to  a  40%  error  (9:  Class  2) 
from  the  actual  aircraft  value.  By  finding  the  most  critical  stabihty  derivatives,  the 
required  level  of  accuracy  in  the  estimation  of  the  stabihty  derivatives  was  determined. 

Straight  and  Level  Fhght. 

The  first  longitudinal  trim  condition  investigated  was  straight  and  level  fhght. 

The  Matlab®  routine  developed  in  Chapter  4  (see  Appendix  B)  was  used  to  determine  the 
sensitivity  of  each  output  to  derivative  change.  Changes  in  control  surface  deflection 
greater  than  ±1.0°  when  the  original  value  was  varied  ±40%  were  noted. 

The  Mowing  stabihty  derivatives:  C^^  ,  C^^ ,  C^^ ,  C^^ ,  C^, ,  C^^^ ,  C„^ ,  C^^ , 

and  C^  had  neghgible  effects  on  control  surface  deflection  values  when  varied  ±40%. 
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Increasing  did  not  produce  a  change  in  5e  larger  than  1.0°,  but  decreasing 
23.8%  from  the  original  value  resulted  in  a  1.0°  decrease  in  5e.  Varying  by  ±31% 

resulted  in  a  change  of  1.0°  from  the  original  value  of  5e.  Overall,  in  straight  and  level 
flight  errors  in  the  longitudinal  stabUity  derivatives  did  not  produce  a  significant  change 
in  the  control  surface  deflection  values. 

Oimbing  Flight. 

Next,  the  longitudinal  stability  derivative  sensitivities  were  investigated  in 
trimmed  climbing  flight.  The  same  stabUity  derivatives  were  evaluated  with  the  same 
criteria  as  straight  and  level  longitudinal  flight.  Changes  in  aU  the  longitudinal  stability 
derivatives  except  had  a  neghgible  effect  on  the  required  control  surface  deflections. 

As  in  straight  and  level  flight,  a  ±31%  error  in  produced  a  1.0°  change  in  5e. 

Overall,  the  errors  in  the  longitudinal  stabihty  derivatives  did  not  produce  large 
errors  in  the  control  surface  deflection  values.  The  accuracy  of  the  longitudinal  values 
was  also  well  known  (reference  10).  Therefore  the  longitudinal  predictions  were  very 
accurate  and  should  match  the  flight  test  results.  In  Chapter  6  the  accuracy  of  the 
longitudinal  predictions  compare  very  well  to  the  measured  flight  test  data. 

Lateral-Directional  Stability  Derivative  Sensitivity. 

The  lateral- directional  stabihty  derivatives  sensitivities  were  found  for  two  cases, 
crabbing  and  steady  turn  flight.  Each  stabihty  derivative  was  varied  ±40%  from  the 
original  value  to  estimate  the  effect  of  errors  in  the  stabihty  derivatives  on  the  predicted 
control  surface  deflections. 
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Crabbing  Flight. 


The  first  lateral- directional  trim  condition  investigated  was  flight  in  a  crabbing 
condition  with  a  constant  sideshp  angle.  Each  lateral- directional  stabihty  derivative  was 
again  varied  ±40%  from  the  original  values.  Changes  of  ±40%  in  the  following  stability 

derivatives  changed  the  control  surface  deflections  less  than  one  degree:  Cy  ,  C^,  ,  , 


,c,,c„  ,C„^,andC„^. 


Small  dihedral  stabihty  derivative,  Cj  ,  errors  produced  large  changes  in  the 


predicted  control  surface  deflection  values.  A  change  of  17.1%  of  C,  in  either  direction 
from  the  original  value  created  a  1.0°  change  in  5  a,  while  a  variation  in  C,  did  not 


change  5r.  Changes  in  the  aileron  control  power  derivative,  C,^  ,  produced  an 


insignificant  change  in  5r,  but  a  18.7%  decrease  from  the  original  value  of  Cj^  or  a  30% 
increase  in  Cj^  changed  5a  1 .0°  from  the  original  value.  Similarly,  varying  the  roU 
couphng  stabihty  derivative,  C,^  ,  had  a  neghgible  effect  on  the  5 r  predictions.  However, 
changing  C,  ±28.6%  resulted  in  a  1.0°  change  in  5a.  Smah  errors  in  the  weathercock 


stabihty  derivative,  C  ,  had  large  ramifications  on  the  predicted  lateral- directional 


control  surface  deflections.  A  ±9.3%  change  in  C^  resulted  in  a  1.0°  change  in  5r.  Also, 
a  ±28.6%  change  in  C^  produced  the  same  error  in  5a.  Finahy,  errors  in  the  mdder 
control  power  derivative,  C^  ,  produced  significant  errors  in  the  lateral- directional 
control  surface  deflection  values.  While  an  increase  in  C„  did  not  effect  5a,  a  22.1% 

■lar 
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decrease  in  from  the  original  value  changed  5  a  by  1.0°.  Rudder  surface  deflection 

was  even  more  sensitive  to  changes  in  C„  .  A  10.4%  increase  and  an  8.6%  decrease  in 

from  the  original  value  produced  a  1.0°  change  in  5r. 

Overall  the  lateral- directional  stabihty  derivatives  were  more  sensitive  to  errors  in 
crabbing  flight  than  the  longitudinal  stabihty  derivatives  in  any  flight  condition.  Errors  in 
lateral- directional  stability  derivative  values  (reference  10),  also  contributed  to  the  errors 
in  the  lateral- directional  axes  predictions.  Therefore  the  open- loop  equations  wiU  not 
predict  required  control  surface  deflection  for  crabbing  flight  as  accurately. 

Steady  Turn  Fhght. 

The  other  lateral- directional  trim  condition  investigated  was  flight  in  a  steady  turn 
with  a  constant  bank  angle.  Each  lateral- directional  stabihty  derivative  was  again  varied 
±40%  from  the  original  values.  Changes  of  ±40%  in  the  fohowing  stabihty  derivatives 
had  a  neghgible  effect  on  the  control  surface  deflection  values:  C„  ,  ,  C„  ,  C„  , 

yp/ryaayar 

C,^,C,,C„^,andC„^. 

As  in  the  crabbing  fhght  case,  the  same  stabihty  derivatives  were  most  sensitive 
to  errors  in  a  steady  turn.  A  decrease  of  1 1 .0%  in  the  original  value  of  C  ^  resulted  in  an 

y  a 

increase  in  5a  of  1.0°.  Eikewise,  only  a  4.8%  decrease  or  a  5.6%  increase  in 

y  a 

produced  an  unacceptable  error  in  5r.  While  errors  in  C,  did  not  effect  5r,  an  error  of 

±6.9%  resulted  in  a  1.0°  change  in  5a.  The  predicted  5r  was  not  influenced  by  errors  in 
either  C,  or  Cj  .  Aileron  deflection  was  significantly  impacted  by  errors  in  the  roh 
coefficient  control  surface  derivatives  though.  An  11.1%  decrease  or  a  14.6%  increase  in 
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resulted  in  a  significant  change  in  5a.  Also,  a  ±15.6%  change  in  produced  a 
similar  error  in  5a.  Small  weathercock  stabihty  derivative,  ,  errors  produced  large 
errors  in  the  control  surface  deflection  predictions.  A  ±5.0%  variation  in  changed  5r 

by  a  degree,  and  if  the  error  increased  to  ±15.4%,  5a  exceeded  the  error  limit.  Finally, 
only  a  ±4.6%  error  in  C„  resulted  in  a  change  of  1.0°  in  5r.  As  the  error  in  C„  was 

decreased  to  - 13.2%  or  increased  to  18.6%  of  the  original  value  of  C„  ,  the  same  error 

^  ^ar 

was  seen  in  5  a. 

As  seen  in  the  crabbing  flight  case,  in  a  steady  turn  errors  in  the  stabihty 
derivatives  produced  unacceptably  large  errors  in  the  lateral- directional  control  surface 
deflection  predictions.  The  higher  uncertainty  in  the  lateral- directional  derivative  values 
also  compounded  higher  inaccuracy  in  the  predictions  when  compared  to  the  longitudinal 
cases.  The  higher  errors  predicted  for  the  lateral- directional  cases  was  observed  in  flight 
test  (see  Chapter  6). 

Summary 

The  stabihty  derivatives  were  prioritized  from  the  most  to  the  least  sensitive  to 
change  in  control  surface  deflection  values.  The  results  are  hsted  in  Table  2  below.  The 
percentage  change  that  produced  a  1.0°  error  in  control  surface  deflection  dictates  the 
level  of  accuracy  required.  The  stabihty  derivatives  were  estimated  in  the  flight  test 
program  to  improve  the  prediction  accuracy  (reference  10).  The  derivatives  in  Table  2 
were  the  focus  of  the  aeromodeUing  tests  to  attempt  to  achieve  the  required  level  of 
accuracy.  In  flight  test,  aU  longitudinal  stabihty  derivatives  were  within  acceptable 
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tolerances.  The  accuracy  of  the  new  and  estimates  did  not  reach  the  accuracy 

required  for  the  lateral- directional  stabihty  derivatives.  Therefore,  the  lateral- directional 
control  surface  deflection  predictions  cannot  be  guaranteed  to  be  within  ±1.0°  of  the 
actual  value.  This  discrepancy  was  observed  in  flight  test  (see  Chapter  6). 

The  prioritized  sensitivities  of  the  stabihty  derivatives  foUow.  A  “■+■”  percentage 
is  an  increase  from  the  original  value,  while  a  is  a  percent  decrease. 
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Table  2.  Stability  Derivative  Sensitivity  Analysis 


STABILITY  DERIVATIVE 

5e 

5a 

5r 

CONDITION 

Rudder  Control  Power  ( C  „  ) 

^ar 

>±40% 

-13.2%/ 

-tl8.6% 

±4.6% 

Banked  Fhght 

Side- force  Coefficient  due  to 
Sideslip  (C^  ) 

>±40% 

-11.0% 

-4.8%  / 
-t5.6% 

Banked  Fhght 

Weathercock  Stability  (C„  ) 

>±40% 

±15.4% 

±5.0% 

Banked  Fhght 

Dihedral  Stabihty  (Cj  ) 

>±40% 

±6.9% 

>±40% 

Banked  Fhght 

Rudder  Control  Power  ( C  „  ) 

“ar 

>±40% 

-22.1% 

-8.6%/ 

-1-10.4% 

Crabbing  Fhght 

Weathercock  Stability  (C„  ) 

>±40% 

±28.6% 

±9.3% 

Crabbing  Fhght 

Aileron  Control  Power  ( C ,  ) 

'aa 

>±40% 

-11.1%/ 

-1-14.6% 

>±40% 

Banked  Fhght 

Roll  Moment  Coefficient  due  to 
Rudder  Deflection  ( C ,  ) 

^ar 

>±40% 

±15.6% 

>±40% 

Banked  Fhght 

Dihedral  Stabihty  (Cj  ) 

>±40% 

±17.1% 

>±40% 

Crabbing  Fhght 

Aileron  Control  Power  ( C ,  ) 

^aa 

>±40% 

-18.7%/ 

-t30.0% 

>±40% 

Crabbing  Fhght 

lift-Curve  Slope  (Cl^  ) 

-23.8% 

>±40% 

>±40% 

Straight  and  Level 
Fh^t 

RoU  Moment  Coefficient  due  to 
Rudder  Deflection  ( C ,  ) 

^  ^ar  ^ 

>±40% 

±28.6% 

>±40% 

Crabbing  Fhght 

Pitch  Coefficient  due  to  Angle 
of  Attack  (C^ , ) 

'  m  a  ^ 

±31.0% 

>±40% 

>±40% 

Straight  and  Level  / 
Climbing  Fhght 
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VI.  Flight  Test  Results 


Overview 

In  Chapter  2,  the  point  mass  model  was  developed  to  provide  a  basis  for  the 
control  laws  for  the  trajectory  of  the  aircraft.  The  outer- loop  equations  were  found  from 
the  nonlinear  equations  of  motion  and  apphed  to  the  actual  Learjet-25  parameters.  In 
Chapter  3,  the  outer- loop  equation  outputs  were  used  with  the  transformation  matrices 
and  the  remainder  of  the  equations  of  motion  to  find  the  control  surface  deflections 
required  to  achieve  the  desired  flight  trajectory  and  orientation.  In  Chapter  4,  the 
equations  developed  in  the  previous  two  chapters  were  used  to  simulate  selected  flight 
trajectories.  By  choosing  the  correct  input  states,  the  required  control  surface  deflections 
for  the  desired  flight  path  were  found.  The  impacts  of  variations  in  the  stabihty 
derivatives  were  investigated  in  Chapter  5  to  determine  which  parameters  require  the 
most  accuracy.  Finally,  in  Chapter  6  the  theoretical  values  tested  using  the  Learjet-25  in 
flight  tests  are  compared  to  the  theoretical  calculations.  The  aeromodel  for  the  Learjet-25 
was  estimated  from  flight  test  to  improve  the  stabihty  derivatives.  Selected  trajectories 
from  the  previous  chapter  were  also  selected  for  verification  in  flight  test.  The  results  of 
the  flight  test  were  compared  to  the  predicted  results  to  determine  the  accuracy  of  the 
model. 

Aeromodel  Data 

Veiidian  oiiginaUy  provided  the  aerodynamic  model  for  the  Leatjet-25.  The 
original  values  were  used  in  the  Matlab®  to  calculate  the  estimated  control  surface 
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deflections.  The  uncertainties  in  the  provided  stabihty  derivatives  could  produce  errors 
larger  than  1.0°  in  the  predictions.  Therefore,  in  the  flight  test  program  the  aeromodel 
parameters  were  estimated  to  produce  improved  values  for  the  stabihty  derivatives 
(reference  10). 

After  stable  data  were  collected  at  each  trim  flight  test  point,  control  surface 
deflection  doublets  were  entered  into  the  aircraft.  Programmed  test  input  was 
automaticahy  input  for  each  of  he,  ha  and  hr  and  the  aircraft  response  recorded.  Using  the 
response  at  each  test  point,  the  stabihty  derivatives  were  estimated  using  Parameter 
Estimation  software  provided  by  TPS.  Average  stabihty  derivative  values  were  found  for 
both  Regime  1  and  Regime  2.  The  effect  of  altitude,  Mach  number  and  eg  location  were 
also  investigated  for  each  stabihty  derivative  (reference  10).  In  general,  the  changes  in 
stabihty  derivatives  were  neghgible  within  each  regime.  The  improved  stabihty 
derivatives  were  then  re- inserted  into  the  Matlab®  routine  and  the  control  surface 
deflections  re-calculated.  The  improved  stabihty  derivative  values  were  used  throughout 
the  study  and  are  listed  in  Table  2.  The  improved  stabihty  derivative  predictions  were 
compared  to  the  flight  test  data  in  the  remainder  of  this  chapter. 

Flight  Test  Data 

Selected  trim  flight  conditions  examined  in  Chapter  4,  were  verified  in  flight  on 
the  Learjet-25  and  compared  to  the  theoretical  predictions.  Trim  conditions  2,  3,  4,  5, 
and  10  were  investigated  in  the  flight  test  program.  Most  of  the  longitudinal  results  (trim 
conditions  1,  2,  and  3)  matched  the  theory  very  closely,  while  lateral- directional  flight 
conditions  (trim  conditions  I,  4,  and  5)  had  more  errors  and  did  not  match  the  predictions 
as  weh.  Each  condition  investigated  during  the  flight  test  program  is  now  examined. 
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The  flight  test  envelope  was  divided  into  three  flight  regimes  to  study  the 
differences  between  low,  middle  and  high-speed  results.  The  three  regimes  were  chosen 
to  investigate  the  effects  of  atmospheric  compressibihty  and  the  non- linear  regions  of  the 
flight  envelope.  The  middle  regime,  Regime  1,  was  defined  as  190  KIAS  to  250  KIAS 
(approximately  400  to  550  ft/s).  Low  speed.  Regime  2,  was  deflned  as  the  stall  speed 
(dependent  on  aircraft  weight,  but  around  120  KIAS)  to  190  KIAS.  The  high-speed 
regime.  Regime  3,  was  from  250  KIAS  to  the  maximum  aircraft  airspeed  (310  KIAS  or 
about  650  ft/s).  Fhght  test  in  high-speed  regime  was  not  completed  due  to  time 
constraints  in  the  test  program. 

The  Learjet-25  had  a  variable  stabilator  (moving  horizontal  tail)  that  changed  the 
flying  characteristics  of  the  airplane  as  the  stabilator  position  was  changed.  Stabilator 
position  was  a  function  of  airspeed  to  maintain  trimmed  flight  Within  each  regime,  the 
stabilator  position  was  held  constant  but  each  regime  had  a  different  trimmed  stabilator 
position.  Therefore,  the  trimmed  aircraft  variables  and  stabihty  derivatives  changed 
between  regimes.  Exact  comparisons  between  regimes  are  not  vahd. 

Straight  and  Level  Flight  Test. 

The  majority  of  the  flight  test  program  was  used  to  examine  trim  condition  I,  and 
aU  the  associated  flight  effects  in  straight  and  level  flight.  Four  flight  conditions  were 
examined  within  trim  condition  1:  altitude  effect,  velocity  effect,  weight  effect  and 
asymmetrical  thmst  effect.  Each  flight  condition  will  now  be  detailed. 
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Altitude  Effects. 


First  the  effect  of  altitude  change  was  investigated  by  selecting  straight  and  level 
test  points  between  9,000  and  20,000  feet  mean  sea  level  (MSL).  The  effect  of  weight 
and  velocity  were  isolated  to  the  greatest  extent  possible  by  flying  each  test  point  with  the 
same  aircraft  configuration  at  each  altitude  in  the  test  band.  The  results  are  shown  in 
Figure  30  through  Figure  33. 

The  altitude  effect  on  the  flight  test  parameters  was  compared  to  the  predictions. 
The  results  were  very  scattered  with  little  correlation.  Therefore  the  data  were  further 
separated  into  three  aircraft  velocity  bands  within  each  fiight  regime.  The  flight  test 
results  in  Regime  1  generally  matched  the  predicted  values,  while  the  Regime  2  flight  test 
data  were  biased  from  the  predicted  curves. 


Figure  30.  Regime  1  Elevator  Deflection  (6e)  vs.  Altitude  (h)  Effect 
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Figure  30  shows  the  he  test  data  in  Regime  1  with  an  overlay  of  the  predicted  he 
for  three  true  airspeed  values.  The  low  speed  region  (V  =  460  ft/s)  described  a  much 
more  shallow  curve  than  theory  predicted.  The  middle  speed  region  (V  =  475  ft/s) 
matched  the  prediction  with  a  +0.2°  bias  for  all  altitudes.  The  high  speed  (V  =  490  ft/s) 
test  points  matched  the  predicted  curve,  with  about  the  same  +0.2°  bias,  but  was  more 
scattered  than  the  mid -speed  data.  Overall  the  flight  test  data  matched  the  predicted  data 
with  a  small  amount  of  variation. 


Figure  31.  Regime  2  Elevator  Deflection  (5e)  vs.  Altitude  (h)  Effect 


In  Regime  2  the  flight  test  data  differed  from  the  predicted  curve  by  almost  +1.0° 
on  average  for  all  three  velocity  bands  (see  Eigure  31).  The  flight  test  data  were  also 
more  scattered  than  in  Regime  1.  Overall  the  data  in  Regime  2  shows  the  same  general 
trend  as  the  predicted  curves,  but  the  correlation  was  not  strong.  The  altitude  effect  on  Ue 
was  very  small,  and  was  masked  by  other  factors,  especially  small  changes  in  airspeed. 
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The  ha  and  hr  were  predicted  to  be  zero  in  trimmed  flight.  Small  values  of  ha  and 
hr  were  measured  in  flight.  These  deviations  result  from  the  impossibUity  to  be  perfectly 
trimmed  with  no  lateral- directional  inputs  during  actual  flight  test.  There  were  no 
significant  values  of  ha  and  hr  in  any  of  the  longitudinal  axes  investigations  data  points. 


Figure  32.  Regime  1  Angle  of  Attack  (a)  and  Pitch  Angle  (0)  vs.  Altitude  (h) 

Effect 


The  altitude  effect  on  a  and  0  was  also  not  very  strong.  Figure  32  show  the  plot 
of  flight  test  data  in  Regime  1  for  a  and  0  respectively,  with  the  predictions  shown  for  the 
same  three  hue  airspeed  bands.  The  data  in  all  three  bands  were  characterized  by  a  lot  of 
scatter  in  the  flight  test  data.  All  data  were  within  2°  of  the  predicted  curves.  The  0 
values  were  closer  to  the  predicted  curves  and  did  not  have  as  much  scatter  as  a. 


81 


Regime  2  Altitude  Effects 
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Figure  33.  Regime  2  Angle  of  Attack  (a)  and  Pitch  Angle  (0)  vs.  Altitude  (h)  Effect 


In  Regime  2,  the  data  were  again  scattered,  increasing  in  error  as  the  altitude 
increased  (see  Figures  33).  The  flight  test  data  were  also  biased  -2°  to  -3°  from  the 
predicted  value  curves.  Also  the  a  values  were  more  scattered  than  the  0  values  at  each 
altitude. 

Overall,  the  altitude  effect  on  the  elevator  surface  deflections  was  very  small. 

The  data  points  were  scattered  by  other  flight  parameter  effects  and  the  predicted  trends 
were  difficult  to  match  to  the  model  predictions.  The  effect  of  altitude  variations  does 
not  have  a  signiflcant  impact  on  the  resultant  errors. 
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In  general,  more  thrust  was  measured  compared  to  the  predicted  curve  (see  Figure 
34).  The  flight  test  data  follows  the  same  general  trend  as  the  predicted  curve.  The 
majority  of  the  flight  test  data  are  within  100  lbs  of  the  predicted  curve.  The  bias  in  the 
data  is  likely  due  to  the  simphfled  thmst  model  used  to  calculate  the  flight  test  thmst. 

The  thrust  was  calculated  from  the  engine  pressure  ratio  and  the  known  area  of  the 
nozzle.  A  more  robust  aircraft  thmst  measurement  should  improve  the  results.  Also  the 
open- loop  equations  assume  the  thmst  is  two  point  forces  acting  at  the  estimated  center  of 
the  nozzle  exit  plane.  This  simplification  also  likely  increases  the  error  in  the  thmst 
prediction. 

Weight  Effects. 

Next  the  effect  of  changing  the  aircraft  weight  and  eg  was  compared  to  the 
predictions.  Since  eg  was  a  linear  function  of  the  Learjet-25  weight  for  aU  test  points,  the 
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control  surface  deflections  were  plotted  as  a  function  of  weight.  Figure  35  shows  the 
combined  results  for  Regime  1  and  Regime  2. 


Figure  35.  Elevator  Deflection  (6e)  vs.  Weight  (W)  Effect 


Eor  Regime  1  the  flight  test  results  were  almost  all  within  ±0.5°  of  the  theoretical 
prediction.  However,  the  flight  test  results  did  not  show  a  signrflcant  relationship 
between  tie  and  variations  in  weight;  actual  ^  was  essentially  constant.  According  to 
theory,  a  3,000- pound  increase  in  weight  would  produce  a  1.0°  decrease  in  ^  in  Regime 
1,  but  this  was  not  observed  during  flight  test. 

In  Regime  2  the  data  were  more  scattered  and  also  did  not  match  the  predicted 
trend.  The  Regime  2  flight  test  results  were  on  average  about  1.0°  higher  than  the 
predicted  value.  Again  the  flight  test  ^  were  nearly  constant  for  all  weights  in  the  test 
envelope.  Theoretically,  a  1,500- pound  increase  in  weight  would  produce  a  1.0°  decrease 


84 


in  ^  in  Regime  2.  Overall,  the  impact  of  aircraft  weight  did  not  have  a  measurable  effect 
on  the  ^  value  from  flight  test  data.  Theoretically  the  impact  was  expected  to  be  small, 
but  not  as  small  as  the  actual  test  data  indicated.  The  significant  difference  between 
Regime  1  and  Regime  2  ^  were  the  result  of  a  change  in  trimmed  stabilator  position. 
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Figure  36.  Angle  of  Attack  (a)  and  Pitch  Angle  (0)  vs.  Weight  (W)  Effect 


The  plot  of  a  and  0  verses  weight  is  shown  in  Figure  36.  Both  Regime  1  and 
Regime  2  are  shown  on  the  plot.  For  Regime  1,  the  flight  test  data  are  scattered  within  2° 
of  the  predicted  curve,  with  all  the  0  data  points  within  1°  of  the  expected  values.  The  a 
values  were  directly  affected  by  the  airspeed  at  which  the  test  point  was  flown.  A  10  ft/s 
change  in  tme  airspeed  resulted  in  an  average  0.5°  to  1.5°  change  in  a.  The  a  and  0 
values  were  more  constant  than  the  predicted  curve.  The  accuracy  of  the  in-flight  a  and 
0  measurements  contributed  to  the  scatter  in  the  data.  The  a  measurement  varied  ±  1 .0° 
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to  2°,  and  the  0  measurement  varied  up  to  ±1.0°.  Therefore,  the  weight  effect  on  these 
flight  parameters  was  not  as  large  as  expected.  The  accuracy  of  the  a  measurement  may 
also  have  been  less  than  for  0.  In  the  aeromodehng  investigation  a  lag  was  found  in  the  a 
and  (3  measurements  (reference  10). 

In  Regime  1  the  flight  test  data  were  consistently  1°  to  2°  lower  than  the  predicted 
values.  The  flight  test  data  showed  the  same  scatter  in  Regime  2  as  in  Regime  1,  with  0 
nearly  hnear  around  the  average  value  and  a  scattered  more  about  the  average  value. 

The  larger  discrepancy  for  a  and  0  in  Regime  2  may  be  due  to  the  fact  that  the  equations 
were  developed  for  Regime  1,  and  small  aerodynamic  non-hnearities  and  the  changed 
position  of  the  stabilator  could  cause  an  error  in  the  predictions. 

The  method  weight  was  measured  is  another  source  of  error  in  the  weight  flight 
test  data.  The  aircraft  weight  was  displayed  in  the  cockpit,  but  was  not  directly  input  into 
the  data  bus  messages  recorded  for  the  flight.  Veiidian  developed  a  routine  that  set  the 
aircraft  weight  at  the  start  of  the  flight  and  decremented  the  weight  at  a  set  rate  according 
to  time  in  flight.  It  was  assumed  that  the  fuel  bum  rate  was  constant,  which  may  not  have 
been  an  accurate  assumption.  Periodically  the  weight  was  updated  from  the  cockpit 
display  to  improve  the  recorded  weight  value.  The  first  few  flights  the  decrement  rate 
was  varied  to  try  to  match  the  aircraft  display  in  the  cockpit.  The  weight  measurement 
error  would  increase  the  scatter  in  the  flight  test  data. 
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Figure  37.  Thrust  (T)  vs.  Weight  (W)  Effect 


Straight  and  Level  Weight  Effects 
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The  weight  verses  thrust  plot  (see  Figure  37)  has  even  more  scatter  than  the 
previous  thmst  plot  (Figures  34).  The  increased  scatter  was  likely  due  to  the  method 
weight  was  measured  (discussed  above),  and  the  same  errors  used  to  calculate  thrust. 
The  data  has  the  same  upward  trend  as  weight  increased,  and  the  predicted  values  are 
lower  than  the  measured  data  as  previously  observed.  Next  the  effect  of  var3dng  the 
velocity  was  investigated. 

Velocity  Effects. 

The  effect  of  changing  the  aircraft  velocity  in  flight  test  is  compared  to  the 
theoretical  prediction.  The  results  for  both  Regime  1  and  Regime  2  are  plotted  in  Figure 
38.  The  flight  test  data  followed  the  same  trend  predicted  by  the  outer- loop  equations, 
with  a  small  bias.  The  data  bias  was  more  pronounced  in  Regime  2. 


87 


Velocity  Effect 
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Figure  38.  Elevator  Deflection  (5e)  vs.  Velocity  (V)  Effect 


In  Regime  1,  the  flight  test  ^  was  consistently  higher  than  the  theoretical 
prediction  by  approximately  0.25°.  AH  flight  test  data  were  within  0.0°  to  0.5°  of  the 
model  prediction.  The  effect  of  velocity  was  not  affected  by  small  changes  in  the  aircraff: 
weight  and  test  point  altitude.  Therefore,  the  effect  of  velocity  on  de  had  the  greatest 
impact  on  the  straight  and  level  flight  results.  In  Regime  2  the  flight  test  de  was 
consistently  0.75°  higher  than  the  predicted  curve.  The  flight  test  data  were  tightly 
grouped  around  the  average  trend- line  and  showed  the  strong  effect  that  a  change  in 
velocity  created  in  the  dg  value.  The  larger  difference  between  predicted  and  flight  test 
data  in  Regime  2  was  similar  to  the  results  as  discussed  above  for  the  altitude  effect. 
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Figure  39.  Angle  of  Attack  (a)  and  Pitch  Angle  (0)  vs.  Velocity  (V)  Effect 


The  a  and  0  results  for  both  regimes  are  plotted  in  Figure  39.  The  flight  test 
results  matched  the  theoretical  predictions  within  ±1.0°,  except  at  the  slowest  airspeeds  in 
Regime  2.  The  data  were  in  a  nearly  continuous  curve  across  both  regimes.  There  was  a 
shght  jump  in  the  theoretical  curve  between  Regime  1  and  Regime  2  of  less  than  0.2°. 
Again  there  was  more  scatter  in  the  a  flight  test  results  than  in  the  0  values.  The  a  and  0 
flight  test  data  matched  the  model  very  closely.  The  velocity  effect  dominated  and  the 
data  were  scattered  by  the  other  factors. 

Overall,  there  was  a  strong  correlation  between  changes  in  tme  airspeed  and  he. 
The  correlation  was  also  seen  in  a  and  0  with  velocity  changes.  Velocity  measurement 
accuracy  was  one  of  the  key  indicators  of  he  prediction  accuracy. 
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The  thrust  measurements  as  a  function  of  velocity  are  shown  in  Figure  40.  The 
flight  test  data  had  the  same  general  trend  as  the  predicted  curve,  but  again  results  were 
biased  up  about  100  lbs.  The  same  problems  experienced  in  the  altitude  investigation 
apply  to  the  velocity  effect  on  thrust. 

Longitudinal  Maneuvers,  Straight  and  Level  Flight  Summary. 

Overall,  the  effect  of  a  produced  a  larger  change  in  the  measured  Ug  than  the 
change  in  aircraft  weight  (and  corresponding  eg).  The  weight  effect  was  masked  by  the 
changes  in  a  from  the  predicted  values.  Airspeed  also  had  a  large  effect  on  ^ 
measurements  in  flight.  The  scatter  in  the  data,  especially  in  Regime  2,  could  be 
explained  by  the  variances  in  the  a  and  V  at  each  test  point  from  the  ideal  condition  used 
for  the  model  calculation.  When  the  identical  flight  condition  parameters  were  used  in 
the  model  the  predicted  values  closely  matched  the  flight  test  data. 
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The  test  point  set-up  had  a  major  impact  on  variations  between  the  flight  test  data 
and  the  theoretical  predictions.  The  flight  test  velocity  had  the  largest  impact  on  changes 
in  the  he  values.  Small  variations  in  test  point  velocity  produced  a  large  change  in  he 
values.  The  efiect  of  weight  variations  was  less  significant,  and  the  effect  of  changing 
altitude  was  very  small.  The  test  set-up  variations  in  (3  and  flight  path  angle  from  the 
specified  test  conditions  were  generally  less  than  ±1.0°  for  all  test  points  and  did  not 
significantly  impact  the  results.  The  bank  angle  varied  up  to  ±8.0°  from  the  specified  test 
conditions.  Often  there  was  a  small  roll  rate  (±0.5°/sec)  that  was  not  observed  in  flight, 
but  impacted  the  flight  test  data.  The  (|)  measurement  also  had  the  largest  oscillations 
about  the  steady- state  value,  up  to  a  ±1 .0° .  The  combined  errors  in  (|)  values  had  a 
significant  impact  on  the  lateral- directional  results.  Combining  variations  resulted  in  a 
composite  change  in  the  control  surface  deflection  value,  which  accounts  for  some  of  the 
scatter  seen  in  the  flight  test  data.  The  thmst  flight  data  has  the  same  general  trends  as 
predicted  by  the  model.  There  was  a  consistent  upward  bias  in  the  flight  test  data  from 
the  expected  curve  though. 

Asymmetrical  Thmst  Effects. 

The  final  trim  condition  I  flight  condition  examined  was  asymmetric  thmst. 

Thrust  from  one  of  the  two  engines  was  progressively  reduced  and  thmst  on  the  other 
engine  was  progressively  increased  while  keeping  the  wings  level  and  the  airspeed  within 
the  required  tolerance.  Data  were  collected  at  each  asymmetric  thmst  setting.  The  data 
were  plotted  as  the  percentage  of  total  thmst  from  the  right  engine. 
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Figure  41.  Aileron  (5a)  and  Rudder  (5r)  Surface  Deflection  vs.  Asymmetric  Thmst  (T) 

Effect 

Limited  data  were  collected,  especially  in  Regime  2.  The  data  were  very 
scattered,  and  did  not  match  the  predicted  values  very  closely  for  lateral- directional 
control  surface  deflection  values.  Figure  41  shows  the  dr  and  da  values  from  flight  test 
compared  to  the  predicted  theoretical  values  for  both  regimes.  The  theory  predicted  a 
larger  asymmetric  thmst  influence  on  dr  (slope  -0.5°  per  10%  change  in  asymmetric 
thrust,  verses  less  than  0. 1°  per  10%  for  the  averaged  values  of  the  flight  test  data  in 
Regime  1).  The  dr  results  were  about  the  same  in  Regime  2.  The  da  results  showed  the 
opposite  trend  result.  The  theory  predicted  about  a  -0. 1°  change  per  10%  change  in 
asymmetric  thmst,  while  the  averaged  flight  test  data  had  a  slope  of  about  -0.5°  per  10% 
change  in  asymmetric  thmst.  The  data  had  a  lot  of  noise  and  stable  points  were  difficult 
to  achieve  in  asymmetrical  conditions.  The  (|)  measurements  also  had  a  lot  of  noise  and 
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were  accurate  to  within  only  a  degree  or  two  of  the  averaged  value  at  most  test  points. 
These  errors  may  account  for  the  inaccuracies  between  the  flight  test  results  and  the 
theoretical  predictions. 

The  de  values  were  compared  to  the  theoretical  predictions.  Theoretically  the  de 
was  constant  for  all  asymmetric  thmst  test  points.  In  Regime  1,  aU  de  values  were  within 
±0.3°  of  the  theoretical  value.  If  the  four  data  points  with  large  changes  in  (|)  were 
discarded,  the  average  flight  test  data  de  was  -0.03°  from  the  theoretical  de  and  aU  flight 
test  points  were  within  ±1.0°  of  the  predicted  value.  In  Regime  2  the  flight  data  de  value 
was  biased  -1.0°  with  very  tittle  scatter  in  the  data. 


Asymmetric  Thrust  Effect 
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Figure  42.  Sideslip  ((3)  and  Bank  Angle  ((|))  vs.  Asymmetric  Thmst  (T)  Elfect 


The  lateral- directional  control  surface  deflection  errors  also  were  increased 
because  the  test  points  were  not  in  perfect  trim  condition  in  most  cases.  Figure  42  shows 
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the  plot  of  a  and  (|),  respectively,  for  each  test  point.  The  theoretical  values  for  the  (|)  were 
expected  to  be  0.0°  because  the  maneuver  was  designed  for  wings  level  flight.  Fhght  test 
points  were  not  perfectly  trimmed  to  wings  level,  straight  flight,  so  this  plot  is  a  measure 
of  the  error  due  to  non- ideal  flight  test  conditions.  The  errors  in  (|)  were  all  less  than 
±4.0° .  The  offsets  from  trim  would  also  contribute  to  the  errors  in  lateral- directional 
control  surface  deflections.  The  (|)  errors  were  largest  for  the  points  around  50% 
asymmetric  T,  which  correspond  to  the  largest  errors  on  the  control  surface  deflection 
plots  (Figure  42).  Also  the  (|)  errors  tended  to  be  more  negative  as  asymmetric  thrust 
approached  0%,  and  more  positive  as  asymmetric  thrust  approached  100%.  This  was 
also  observed  in  the  da  and  dr  values  at  the  extreme  asymmetrical  thmst  conditions. 

Lateral- Directional  Maneuvers,  Straight  and  Level  Flight  Summary. 

The  lateral- directional  stabihty  derivatives  had  a  much  larger  uncertainty  than  the 
longitudinal  stabihty  derivatives.  This  lateral- directional  stabihty  derivative  uncertainty 
contributed  to  variations  in  model  control  surface  deflections  as  indicated  by  the 
sensitivity  analysis.  The  most  critical  stabihty  derivatives  ,  and  ) 

created  a  large  error  for  a  smah  variation.  Finally,  a  bias  was  observed  in  both  the  da  and 
dr  values.  For  instance,  in  Flight  10  a  constant  -1-1.6°  bias  in  dr  was  observed  in  ah 
straight  and  level  test  cases.  The  da  also  had  a  smah  bias  of  about  0.3°  for  Flight  10.  The 
da  and  dr  were  manually  trimmed,  and  the  smah  error  could  not  be  detected  from  the 
cockpit.  The  trim  error  contributed  directly  to  an  error  in  the  control  surface  deflection 
results. 
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Overall  the  asymmetric  thmst  flight  test  data  did  not  match  the  theoretical 
predictions.  The  correct  trend  was  seen  in  the  data,  but  the  flight  test  results  had  large 
differences  from  the  expected  values.  The  difflculty  in  estabhshing  trimmed  flight  in  an 
asymmetrical  flight  condition  and  the  inaccuracy  in  the  bank  angle  measurement 
aggravated  the  errors.  Next  trim  conditions  2  and  3  were  investigated  jointly  in  the  flight 
test  program. 

Climbing  and  Descending  Flight  Test. 


Figure  43.  Elevator  Deflection  (5e)  vs.  Flight  Path  Angle  (y)  Effect 


The  effect  of  flight  path  angle  changes  was  investigated  for  both  climbing  (trim 
condition  2)  and  descending  (trim  condition  3)  flight.  The  y  effect  on  de  is  shown  in 
Eigure  43  for  both  regimes.  The  y  effect  on  he  was  predicted  to  be  very  small,  and  the 
flight  test  data  matched  the  predictions.  In  Regime  1,  all  the  flight  test  data  were  within 
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±0.25°  of  the  theoretical  curve.  In  Regime  2  the  flight  test  he  was  biased  +0.7°  from  the 
predicted  value.  The  flight  test  data  were  also  more  scattered  about  the  average  value  in 
Regime  2  than  in  Regime  1 .  Overall,  the  predicted  y  correlation  to  he  matched  flight  test 
results,  but  the  impact  of  y  on  he  was  very  small. 


Figure  44.  Angle  of  Attack  (a)  and  Pitch  Angle  (0)  vs.  Flight  Path  Angle  (y)  Effect 


Figure  44  shows  the  comparison  of  flight  test  y  effect  on  a  and  0  respectively  to 
the  theoretically  predicted  relationship.  In  Regime  I,  all  flight  test  data  for  both  a  and  0 
were  within  ±1.0°  of  the  model  prediction  for  all  y  values.  In  Regime  2  the  flight  test 
data  were  within  ±2.0°  of  the  predicted  value,  with  a  slight  bias  of  about  -1.0°.  As  seen 
for  other  longitudinal  investigations,  the  scatter  in  the  a  data  points  was  larger  than  the 
scatter  in  the  0  values.  Overall,  the  flight  test  data  correlated  with  the  predicted  values 
for  aU  y.  The  impact  of  y  was  not  significant,  but  was  predicted  properly  by  theory. 
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Figure  45.  Thrust  (T)  vs.  Flight  Path  Angle  (y)  Effect 

The  flight  path  angle  had  a  significant  impact  on  the  thrust  required  to  maintain 
trimmed  flight.  Figure  45  shows  the  thmst  data  collected  as  flight  path  angle  was  varied. 
The  flight  test  roughly  had  the  same  trend  as  the  predicted  curve,  but  as  data  was  taken 
further  from  the  center  test  condition,  the  difference  between  the  measured  and  predicted 
data  grew  large.  The  flight  test  data  followed  a  nearly  linear  trend,  while  a  trigonometric 
function  was  predicted  from  the  model.  The  increased  error  at  the  extremes  of  the  thrust 
envelope  was  likely  due  to  errors  in  the  calculated  thmst  from  flight  test  and  the  model. 

Crabbing  Flight  Test 

Next  a  lateral- directional  trim  condition  was  investigated,  the  effect  of  a  non- zero 
sideslip  angle  in  level  flight.  With  (3  specified,  the  (|)  required  to  maintain  y=  H  =  0°  was 
determined.  The  test  was  limited  to  ±10°,  due  to  aircraft  safety  concerns. 
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Figure  46.  Aileron  (5a)  and  Rudder  (5r)  Surface  Deflections  vs.  Sideslip  Angle  ((3)  Effect 

Regime  1  and  Regime  2,  da  and  dr  were  compared  to  the  da  and  dr  predicted  by  the 
model.  The  dr  and  da  are  shown  in  Figure  46.  The  flight  test  data  for  da  and  dr  were 
linear,  as  was  predicted  by  the  model.  In  addition,  the  slopes  of  the  flight  test  data  had 
similar  slopes  to  that  predicted  by  the  model.  When  the  slopes  did  vary  from  the  model 
predictions,  the  slope  of  the  flight  test  data  was  higher,  but  only  by  approximately  0.1°  of 
d  per  degree  of  (3 .  It  was  also  noted  that  at  [3  =  0°,  da  was  negative  and  dr  was  positive. 

The  deflection  at  (3  =  0°  had  much  more  effect  on  moving  the  data  away  from  the 
prediction  than  the  slopes  of  the  data.  The  rudder  showed  zero  deflection  at 
approximately  (3  =  -2°,  and  the  ailerons  showed  no  deflection  at  approximately  a  =  -0.5°. 

It  is  unknown  if  this  was  an  instrumentation  issue,  a  shghtiy  bent  aircraft,  non- zeroed 
trim,  or  a  combination  of  each  factor. 
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Sideslip  Angle  Effect 
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Figure  47.  Bank  Angle  ((|))  vs.  Sideslip  Angle  ((3)  Effeet 


The  (|)  flight  test  data  were  also  graphed  for  the  steady- heading  sideshps,  as  seen 
in  Figure  47.  For  Regime  1  and  Regime  2,  the  slope  from  the  flight  test  data  was 
approximately  2.5  times  greater  than  the  predieted  slope. 

Overall,  the  data  quality  from  flight  test  was  good  and  showed  a  hnear 
relationship  between  eontrol  surface  deflection  and  |3,  and  (|)  to  |3.  The  slopes  of  the 
control  surface  deflection  to  a  closely  matched  model  prediction,  while  the  slope  of  (|)  to  a 
differed  from  the  model  by  a  factor  of  2.5. 

Level,  Steady  Turn  Flight  Test. 

The  flight  test  results  for  trim  condition  5,  did  not  match  the  predictions.  The 
open- loop  equations  predict  a  certain  path  with  either  a  specifled  bank  angle  or  sideshp 
angle.  When  (|)v  was  specifled,  the  predicted  bank  angle  was  over  the  aircraft  limits, 
which  required  a  much  larger  da  and  dr  than  measured  in  flight  test.  The  aircraft  was 
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flown  with  nearly  zero  sideslip,  but  when  (3=0°  was  specified  the  open- loop  equations 
predicted  very  large  control  surface  deflection  to  command  the  turn,  which  limited  the 
maximum  (|)v.  In  flight,  bank  angle  is  used  to  command  a  turn,  which  requires  little 
control  surface  deflection.  The  flights  should  be  repeated  attempting  to  recreate  the  exact 
flight  conditions.  The  flight  trajectory  flown,  did  not  match  the  mathematically  predicted 
trajectory.  Also  turn  rate  had  a  major  impact  on  the  predicted  control  surface  deflection 
values.  The  turn  rate  was  not  specified  for  the  flight  tests,  and  was  not  measured.  The 
predicted  turn  rate  should  be  compared  to  the  actual  flight  test  turn  rate  to  make  a  valid 
comparison. 

Pitch-Over  Flight  Test. 

Trim  condition  1^  flight  trajectory  also  did  not  match  the  predicted  results.  The 
predicted  trajectory  was  for  straight  and  level  flight,  with  a  pitch  angle  rotation  about  the 
yb-axis  created  by  deflecting  de  only  in  Ig  flight.  To  achieve  the  same  pitch  rotation  rate 
in  flight  test,  the  maneuver  was  entered  in  a  climb  with  a  pitch- over,  and  less  than  Ig 
flight.  The  lower  the  load  factor,  the  higher  the  pitch  rate.  Post- flight,  load  factor  and 
pitch  rate  were  varied  to  attempt  to  re-create  the  flight  trajectory.  The  equations  did  not 
predict  the  same  trajectory  flown,  so  the  required  he’s  did  not  match. 

Summary 

The  longitudinal  flight  test  data  provided  a  good  vahdation  of  all  the  model 
predictions,  except  trim  condition  10.  In  straight  and  level  flight,  velocity  was  the 
dominant  input  parameter  for  determining  the  required  de.  Weight  had  a  small  effect  and 
altitude  an  even  smaller  impact.  The  flight  path  angle  had  a  minimal  impact  on  the  de. 
Thmst  setting  had  the  greatest  impact  on  the  flight  path  angle,  de  was  not  significantly 
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changed  by  changes  in  y.  Trim  condition  W  could  not  be  reproduced  as  a  trim  condition 
using  the  open- loop  model. 

The  lateral- directional  flight  test  data  provided  limited  vahdation  of  the  model 
predictions.  The  da  and  dr  curves  had  the  same  general  trends  for  both  the  asymmetrical 
thmst  and  sideshp  investigations.  The  data  was  biased  by  errors  in  the  lateral- directional 
stabihty  derivatives  and  biases  in  the  lateral- directional  measurements  on  the  aircraft. 

The  lateral- directional  data  could  not  be  guaranteed  to  be  accurate  to  within  1.0°.  The 
trim  condition  5  flight  trajectory  could  not  be  reproduced  using  the  model  predicted 
trajectory. 
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VII.  Conclusions  and  Recommendations 


Conclusions 

Model  Development. 

Outer- loop  equations  were  developed  to  produce  control  laws  for  an  aircraft. 
Nested  loop  equations  were  developed  for  the  model.  In  the  outer- loop  controller,  the 
aircraft  was  modeled  as  a  point  mass.  All  aircraft  forces  were  balanced  to  create 
equations  to  solve  the  control  variables.  For  specified  state  variables  and  either  bank 
angle  or  sideslip  angle  to  define  the  flight  trajectory  the  control  variables  were 
determined. 

The  inner- loop  uses  the  control  variables  from  the  outer- loop  to  calculate  the 
required  control  surface  deflections.  In  the  inner- loop,  the  six-degree  of  freedom  model 
was  used  to  balance  the  moments  on  the  aircraft.  With  the  control  variables  determined 
from  the  flight  trajectory  in  the  outer- loop,  the  inner- loop  calculates  the  required  5a,  6e, 
and  5r  to  achieve  the  desired  aircraft  orientation. 

Trim  Conditions. 

Ten  trim  conditions  were  then  defined  to  allow  simphfying  assumptions  to  the 
model  equations.  Trim  condition  I,  is  straight  and  level  flight.  Trim  condition  2  is 
defined  as  a  straight  and  level  climb.  A  straight  and  level  descent  describes  trim 
condition  3.  A  non- zero  sideslip  angle  in  level  flight  is  trim  condition  4.  Trim  condition 
5  is  a  steady,  level  turn.  By  adding  a  positive  flight  path  angle,  trim  condition  6  describes 
a  steady  climbing  turn.  Similarly,  trim  condition  7  is  a  descending  turn.  Trim  condition 
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8  is  a  barrel  roU.  Next,  for  trim  condition  9,  an  AC- 130  gunship  orbit  is  developed. 
Finally,  in  trim  condition  1^,  a  steady  pitch- over  maneuver  is  described.  Model 
equations  were  developed  for  each  trim  condition. 

Model  Implementation. 

The  equations  were  then  implemented  using  a  Learjet-25  aircraft  aeromodel  in  a 
Matiab®  routine.  Each  trim  condition  was  plotted  to  determine  the  relationships  between 
the  relevant  input  variables  and  the  outputs.  The  impact  of  each  input  variable  on  the 
output  for  the  Learjet-25  was  determined.  The  trim  condition  graphs  were  used  to 
determine  which  input  variables  were  focused  on  during  the  flight  test  investigation. 

Predicted  Results. 

Velocity  had  the  largest  influence  on  trim  condition  I  results  with  weight  and 
altitude  having  small  impacts  on  control  variables  and  5e.  For  asymmetrical  thrust,  there 
were  linear  relationships  for  both  5  a  and  5  r  as  the  difference  in  thmst  between  engines 
was  increased.  For  trim  condition  2  and  3,  thrust  was  the  most  important  parameter.  The 
flight  path  angle  was  limited  by  the  thrust  for  both  climbing  and  descending  flight,  hr 
trim  condition  4,  a  hnear  relationship  was  again  found  for  both  5  a  and  5  r  as  the  sideslip 
angle  was  increased.  For  trim  condition  5,  the  steady  turn  was  defined  as  a  constant  turn 
rate  with  a  load  factor  of  1.0  to  allow  a  vahd  comparison.  To  simulate  an  actual  turn,  the 
load  factor  should  be  varied  with  bank  angle,  and  (3  should  be  held  to  0.0° .  Again  a  linear 
relationship  was  observed  for  both  5  a  and  5r  as  the  bank  angle  was  increased.  For  trim 
conditions  6  and  7  5a  and  5r  had  relationships  nearly  identical  to  trim  condition  5.  The 
only  difference  between  the  two  trim  conditions  was  the  required  thmst  was  dependent 
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on  the  flight  path  angle  similar  to  trim  condition  2  and  3 .  For  the  barrel  roll,  a  constant 
roU  rate  was  chosen,  and  (3  was  assumed  to  be  zero.  This  flight  condition  describes  an 
aileron  roU.  To  offset  from  the  roU  axis,  the  sideshp  angle  must  be  increased.  For  trim 
condition  8,5a  was  the  limiting  factor.  Trim  condition  9  was  a  subset  of  trim  condition  5. 
Formulas  were  found  relating  height  above  the  target,  weapon  depression  angle,  and  orbit 
radius.  Finally  a  pitch- over  maneuver  was  examined.  The  pitch  rate  was  set,  and  load 
factor  varied  from  1.0  to  0.4  g’s.  The  pitch-over  was  limited  by  the  maximum  elevator 
deflection. 

Predicted  Derivative  Sensitivity. 

The  equation  predictions  are  very  sensitive  to  errors  in  the  stabihty  derivatives. 
The  most  sensitive  stabihty  derivatives  are  C„  ,  C„  ,  C„  ,andC,  .  An  error  greater 

than  seven  percent  in  any  of  these  stabihty  derivatives  creates  a  1.0°  difference  in  the 
predicted  control  surface  value.  To  get  the  desired  accuracy  in  the  equation  predictions, 
the  stabihty  derivatives  must  be  known  to  great  certainty.  The  stabihty  derivatives 
originahy  provided  by  Veiidian  had  as  large  as  40%  errors  from  the  values  found  in  flight 
test.  The  longitudinal  stabihty  derivatives  had  been  determined  in  previous  studies  and 
were  much  more  accurate  than  the  lateral- directional  stabihty  derivatives  (reference  10). 
Unfortunately  the  lateral- directional  stabihty  derivatives  required  the  most  accuracy 
according  to  the  sensitivity  analysis.  Cj^  and  Cj  were  not  determined  to  within  the 

required  accuracy.  Therefore  the  lateral- directional  predictions  cannot  be  guaranteed  to 
be  within  1.0°  of  the  flight  test  results. 
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Flight  Test  Results. 

The  equation  predictions  were  partially  verified  by  the  fiight  test  program.  The 
longitudinal  test  results  matched  the  predictions  for  trim  conditions  I,  2  and  3 .  The 
lateral- directional  test  results  had  errors  from  the  predictions  for  trim  conditions  !_  and  4, 
but  were  within  acceptable  limits.  Trim  conditions  5  and  10  test  predictions  did  not 
match  the  predictions.  The  errors  were  probably  due  to  not  flying  these  maneuvers  at  Ig 
and  errors  in  the  stabihty  derivatives.  Also  the  last  two  trim  conditions  could  not  be 
exactly  reproduced  by  the  same  flight  trajectory  using  the  control  equations. 

The  thrust  predictions  had  large  errors  and  were  not  very  accurate.  The  thrust 
predictions  had  two  major  sources  of  error.  First,  the  aircraft  measurement  and  model  for 
thrust  introduced  an  error.  Thrust  was  calculated  by  measuring  the  engine  pressure  ratio 
and  the  pressure  at  the  compressor  face.  The  pressure  at  the  nozzle  was  calculated  and 
multiphed  by  the  nozzle  area  to  estimate  the  force.  The  thmst  estimate  was  not  very 
accurate.  From  flight  test,  the  engines  were  balanced  from  all  readings  in  the  cockpit,  but 
post- flight  up  to  a  200  lbs  difference  between  engines  was  noted  at  some  test  points. 

Also,  test  points  flown  at  nearly  identical  test  conditions  on  different  flights  had  up  to  100 
lbs  difference.  The  second  source  of  error  was  that  the  thmst  modeled  as  a  point  force 
acting  at  the  estimated  center  of  the  nozzle.  The  model  simphfications  may  have  also 
introduced  an  error. 

Error  Sources. 

The  equations  were  able  to  accurately  predict  non-aggressive  level  flight 
maneuvers,  which  would  be  the  majority  of  flight  for  a  UAV.  The  complex  maneuvers 
were  not  vahdated  in  the  flight  test  program.  The  steady  turns  predicted  by  the  equations 
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were  also  not  reproduced  in  flight.  Steady  turns  are  a  balance  of  turn  rate,  bank  angle  and 
sideslip  angle.  The  equations  predict  a  different  combination  of  these  three  parameters 
than  were  flown.  The  predicted  steady  turn  may  be  a  vahd  solution  to  the  turn,  but  it  was 
not  verified  in  flight  test.  The  control  surface  deflections  are  smaller  than  predicted  due 
to  the  use  of  bank  angle  to  balance  the  forces.  The  pitch- over  maneuver  also  was  not 
vahdated  in  flight.  Pointing  the  bottom  of  the  aircraft  at  a  stationary  target  on  the  ground 
was  predicted  by  the  equations,  but  may  not  be  possible  on  an  actual  aircraft.  Further 
tests  and  research  is  needed. 

The  model  was  accurate  near  the  center  test  point  condition,  but  the  errors  grew 
larger,  the  farther  from  the  equihbiium  test  point.  The  errors  have  three  primary  sources. 

Stabihty  Derivative  Variabihty. 

First,  the  stabihty  derivatives  were  a  function  of  multiple  variables.  From  flight 
test,  as  Mach  number,  altitude  and  dynamic  pressure  changed  the  stabihty  derivative 
values  changed.  Temperature  and  air  density  also  had  a  small  impact  on  the  predictions. 
To  predict  the  required  control  surface  deflection  with  great  accuracy,  the  stabihty 
derivatives  would  need  to  be  mapped  for  the  entire  flight  envelope  to  the  desired 
accuracy  determined  in  the  sensitivity  analysis. 

Instmmentation  Errors. 

Second,  the  aircraft  measurement  instmmentation  needs  to  be  improved.  Small 
errors  in  the  measured  angle  of  attack,  sideshp,  bank  angle  and  airspeed  directly  impacted 
the  accuracy  of  the  flight  test  results.  It  was  also  noted  that  altitude  measurement 
changed  when  maneuvers  other  than  straight  and  level  flight  were  flown.  The  altitude 
error  was  probably  due  to  no  correction  in  the  static  pressure  measurement  when  flow 
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over  the  static  port  was  disrupted  in  other  than  straight  and  level  flight.  A  Pitot  static 
cahbration  should  be  performed.  The  accelerations  and  angular  rotation  measurements 
were  also  very  noisy,  which  resulted  in  an  averaged  value  with  occasional  large 
uncertainties.  To  achieve  the  accuracy  needed  for  the  model  equations,  very  accurate 
measurements  are  required.  Finally  there  was  no  direct  method  to  record  aircraft  weight. 
The  data  was  displayed  in  the  cockpit,  but  the  data  recorded  from  the  data  bus  was  an 
estimate  provided  by  Veiidian. 

Non- linear  Effects. 

Finally,  at  the  edges  of  the  flight  envelope  non-hnear  flight  effects  become 
evident.  As  stall  is  approached,  the  lift  curve  slope  becomes  non-hnear.  Also  at  higher 
Mach  numbers,  compressibility  effects  must  be  considered.  At  higher  dynamic  pressures 
the  assumption  of  a  rigid  body  may  not  be  valid.  Another  concern  not  addressed  is  the 
movement  of  the  center  of  gravity  as  fuel  was  burned.  These  effects  are  not  included  in 
the  model  and  therefore  would  introduce  an  error.  The  model  accuracy  would  decrease 
the  further  from  the  center  test  condition  the  maneuver  was  flown.  FinaUy  thmst  was 
assumed  to  be  invariant  as  altitude,  velocity  and  temperature  changed.  This  simphfying 
assumption  likely  introduced  an  increasing  error  as  test  data  was  coUected  away  from  the 
center  test  condition. 

Recommendations 

The  open- loop  equations  provide  a  good  prediction  of  the  required  control 
settings  for  some  of  the  specified  trim  conditions.  With  additional  work  the  model  could 
be  implemented  on  a  UAV  apphcation.  This  section  will  address  shortcomings  identified 
during  research  and  additional  issues  that  need  to  be  addressed. 
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Validate  Turn  Predictions. 


The  steady  turn  predictions  need  to  be  validated.  Post- flight  test,  the  effect  of 
load  factor  was  added  to  the  equations  to  predict  a  loaded  turn.  The  equations  stiU 
predicted  much  higher  control  surface  deflections  than  were  seen  in  flight.  Pilot’s 
typically  use  more  bank  angle  than  control  surface  deflection  to  turn  the  aircraft.  The 
equations  predict  a  small  bank  angle  and  large  control  surface  deflection  trajectory.  The 
predicted  trajectory  and  control  surface  deflections  should  be  flight  tested  to  validate  the 
predicted  solution.  The  predicted  trimmed  turn  may  be  a  vahd  turn  solution,  just  not  the 
ones  found  in  flight  test. 

Improve  Thrust  Model. 

The  thmst  model  must  be  improved.  The  model  provided  for  flight  test  did  not 
provide  the  accuracy  to  totally  validate  the  model  thmst  predictions.  Near  the  center  test 
conditions  the  results  were  good,  but  at  the  upper  end  of  the  thmst  envelope  the  flight  test 
measurements  differed  significantly  from  the  predicted  results.  The  thmst  model  had 
many  simplifying  assumptions,  which  may  need  to  be  readdressed  to  accurately  predict 
the  proper  engine  settings. 

Include  Wind  Effect  in  Model. 

The  uncertainties  in  the  atmosphere  need  to  be  included  in  the  model.  Wind  gusts 
and  other  variations  in  wind  should  be  included  to  better  model  the  actual  flight 
conditions.  A  wind  gust  would  add  a  force  and  moment  to  the  system,  which  is  currently 
not  included  in  the  force  and  moment  summations.  The  wind  variations  would  dismpt 
the  assumption  of  zero  net  forces  and  moments,  and  the  aircraft  would  no  longer  be  in  a 
trimmed  flight  condition. 
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Investigate  Transitions  between  Trim  Conditions. 

The  model  assumes  stable,  trimmed  flight  conditions.  No  control  input  is 
provided  to  get  to  the  trimmed  condition.  Typically  large  inputs  are  commanded  to  move 
the  aircraft  to  a  new  flight  condition.  To  transition  from  one  trimmed  condition  to  the 
next,  additional  control  commands  must  be  input.  If  the  trim  condition  control 
commands  are  the  only  inputs  into  the  aircraft,  the  aircraft  may  slowly  approach  the  trim 
condition,  but  may  never  reach  the  desired  trimmed  condition.  The  transition  between 
trimmed  conditions  was  not  addressed  in  the  open- loop  model.  Control  surface  limits 
may  also  be  exceeded  during  transitions  between  trim  conditions. 

Determine  Stability  Derivatives  Accurately. 

The  stabihty  derivatives  require  very  accurate  estimates  to  make  valid  predictions. 
Wind  tunnel  estimates  of  stabihty  derivatives  do  not  produce  the  required  stabihty 
derivative  accuracy.  An  extensive  flight  test  program  is  required  to  provide  an  accurate 
aeromodel  of  the  aircraft. 

To  implement  the  control  loop  equations  on  a  UAV,  very  accurate  flight  condition 
measurements  would  be  required.  The  test  aircraft  had  some  of  the  best  measurement 
instmmentation  economicahy  available.  The  measurements  were  still  near  the  acceptable 
accuracy  limits  and  exceeded  error  limits  in  some  cases.  Extremely  accurate  angle, 
velocity  and  acceleration  measurements  would  be  needed  to  implement  the  open- loop 
control  model  on  an  aircraft. 

Measure  Weight  Accurately. 

The  weight  measurement  was  estimated  since  the  flight  data  was  not  available  on 
the  data  bus.  The  weight  was  only  available  as  a  display  in  the  cockpit.  A  constant  fuel 
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bum  rate  was  assumed  and  aircraft  weight  was  decremented  as  a  function  of  flight  time. 
The  decrement  rate  was  varied  the  first  few  flights  to  find  the  most  accurate  weight 
estimate.  Also  the  weight  was  periodically  reset  to  the  cockpit  display  to  increase 
accuracy.  The  weight  value  occasionally  drifted  from  the  actual  value  though.  The 
weight  of  the  aircraft  should  be  accurately  measured  to  provide  the  accuracy  needed  to 
validate  the  weight  related  model  predictions. 

Automate  Matlab®  Routine. 

The  Matlab®  routine  could  be  automated  to  accept  flight  condition  inputs  and 
output  the  required  control  settings.  Shell  programs  could  be  developed  to  take  inputs 
directly  from  the  air  data  computer  for  air  density,  temperature,  current  velocity,  and 
aircraft  weight.  The  shell  program  could  be  designed  to  prompt  for  the  velocity,  flight 
path  angle,  and  either  bank  angle  or  sideshp  angle  for  the  desired  trajectory.  Additionally 
rotational  rates  and  angular  accelerations  could  be  input  for  the  three  axes.  The  output 
control  surface  deflections  and  engine  setting  could  be  sent  directly  to  the  aircraft  bus 
controller.  Also  the  angle  of  attack,  sideshp,  and  bank  angle  could  be  compared  to  the 
aircraft  measurements  to  monitor  the  aircraft  trajectory. 
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APPENDIX  A:  SAMPLE  CALCULATIONS 


Introduction 

This  appendix  includes  sample  calculations  to  support  the  equations  developed  in 
the  inner-  and  outer- loop  equations.  The  numbers  in  parenthesis  correspond  to  the 
equation  number  in  the  main  text  of  the  document. 


Ligure  48.  Lorce  Definitions 


Outer-Loop  Calculations 

The  outer- loop  equations  were  redefined  in  Chapter  2.  The  derivations  for  the 
equations  follow.  Lirst  the  forces  from  Ligure  48  were  summed  in  all  three  axes 

EL^  =T-D-nmgsinY  =  max  =  m’^ 

T  —  D 

- ng  sin  a  (Eqn.  5) 

m 

L  L^^  =  -L  cos  (|)v  -  Ly  sin  (|)v  -I-  n  m  g  cos  Y  =  m  az  =  m  (- V  5<) 


111 


Lcos4>, +Fsin4>^  ng 

^  l<= - - ^cosa 

mV  V 

(Eqn.  6) 

L  =  Fy  cos  (|)v  -  L  sin  (|)v  =  m  ay  =  m  (V  cos  y  l¥) 

_  -Lsin4>^ +FyCOs4>, 

mVcos  a 

(Eqn.  7) 

The  following  forces  were  defined  in  Chapter  2 

^  ^La 

Vo 

(Eqn.  1) 

oc+Cdo  ) 

(Eqn.  2) 

-  S 

Fy=q  ^  P 

y  ^  v"  s 

(Eqn.  3) 

T=mgp 

(Eqn.  4) 

The  defined  forces  were  substituted  into  equations  (5)  through  (7) 
a+C  )-ngsin  a 

mV„ 


,  cos(|)C,  (a+a  )h — (3 sin 6 

— — - - ^cos  a 

V  V 


.  — Bcostl),  -  sintl)  C,  (a+a  ) 

V  cosy 

Simphfying  assumptions  were  then  made  for  each  trim  condition. 
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Trim  Conditions  1,  2,  3,  and  4. 

Trim  condition:  S<=  0  °  ((3  independent) 


0  =  ^  Bcos4>,-sin4>C  (a+a„  ) 


[Sw 


(from  Inequation) 


c 

— Cy  |3 

^  =  .an.^.=  ^  ' 


cost])^  Cl  (a+a„ 


tan4>,CL  (a+a„  ) 


c  yp 

^\\r 


_ ^  2  /  S  ^ 

0  =  coscP^Cl^^  (oc+O;,^  )h — —Cy  |3sin(l)^  -rung  cos  a  (from '^■equation) 


c  yp 

Ow 


nmgcosa  .  ^  (  )  S 

— —  =  cos6C,  (a+a  j+— ^C„  |3sin6 

q-y2  L„  V  S  yp  ”  ' 


(substituting  for  P) 


nmgcosa  ^  t^ii^vCL^(oc+a  ) 

— —  =  cos6C,  la+a„  In - C„  - — -^ - ^sin6 


rung  cos  a  ^  ^  sin  ^  (t), 

=C,  io^+o^o,  J  cos(^v+ - r 


nm^osa^^^^v  =  Cl^  (a+ )(cos>,  +  sin  =  Cl^  (oc+a„^ ) 


1  X^nmg  cos  a 


costp^  -a„ 
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—\T  2 

0  =  gf^--^-7(CD„oc+C  )-ngsin  a 
mV„ 


(from  ’^equation) 


u=  ^  ^  (Cn  a+Cn  )  +  nsina 
mgV„"  ^  “  ° 


(substituting  for  a) 


yc„„^+C„  +n.ma 

mgV„^  C, 


C  — 

M'=7;^(ncosacos(|>J+-3-— (Cp  -Cq  a  )+nsina 


mgV„ 


Returning  to  the  original  equation,  dividing  by  ^  ,  and  simplifying 


4>vC  L  Joe + a„  J  = C  [3  cos 

Ow 


^C^^|3cos4>, 


CL„loc+a„^, 


(substituting  for  a) 


sin  (])„=■ 


^C^^|3cos4>, 

_ _ 

1  fV^^nmgeosa 
Cr  - — ; - COS4>, 


-Oto,+OCo, 


sin(|)v  = 


C,,S.q 

V^^nmgS  „  cos  a 


(Eqn.  8) 


From  the  trigonometric  relationship 


— C,^|3qV^ 

cos^t])^  =l-sin^(l)^  =  1-  ^ - 

nmg  cos  a 
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COS(^v=.  1-  ^ 


C^.S.qV^a 


nmgS  „  cos  a 


(Eqn.  9) 


Substitute  cos(|)v  into  a 


1  V„nmgcosa  ,  '^y 

a  = - 2 — -  1-  __ 


Cy  S.qV^a 


nmgS  „  cos  a 


V^^nmgcosa^ 


(Eqn.  10) 


Substitute  cos(|)v  into  |i 


-  ,  f  C.^S.qV^a  ^  ^  . 

ncosa^l-  ^ ^  -Co^J+nsina 


V^nmgS^cosa  mgV^ 


U=— i - -ICn  -Cn  a„  1  +  ^  - ^ 

mgV„"  “  D.  o,;  . 


Co.ncosa]  Cj^C^.S.qV^a 


C,  V„WgS 


+  n  sin  a 


(Eqn.  11) 


Trim  condition  ]_:  S<=  l^F^  0  °  ((|)  v  independent)  (from  Equation  8) 


C,S.q 


VjnmgS^  cos  a 


C,S,q  V* 


(Eqn.  12) 


Erom  previous  equations 
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(Eqn.  13) 


(  \T^ 


a  = 


C, 


nmg  cos  a 


cosc])^ 


oc 


C. 


{a=  cos  a  cos  ([)',)+  ^  ,  (Cp  -Cq  )+nsina  (Eqn.  14) 


C, 


mgV' 


Trim  Conditions  5,  6,  7,  and  9. 

Trim  condition:  5<=  0  °  and  1^  =  05>j,  (or  H(t)  =  05^  t) 


—\T  2 

0  =  gfa--^^(c  a+C  )-ngsin  a 
mV„ 


(from  ’^equation) 


0^ 


cos4>,Cl^  (a+a„^ )+  |3sin 


'w 


-  nmg  cos  a  (from  '^■equation) 


mV 03q,  cosy  = 


|3cos4>,  -  sin  (oc+a„^ ) 


(from  Inequation) 


Erom  the  second  and  third  equations 


cos 


4>vC,  (a+a„  )  +  ^C^  Bsin(^,  =- 


V; 


«  ^yp  Yv  “ 


[3  cos4>,  -  sin  {a+ )  =  ^^mV05>p  cosy 


c  yp 

Ow 


qV^ 


Cl  (a  +a„  )  — ^C  p 

*-cc  v  ol  g  yp  ^ 


'w 


-C  P  -CL(a+a  ) 


cos  d, 

Sind, 


-mcos  a 


ng 

VCl5„ 


sint]). 


Vf 


CL/(a+a„ 


+ 


—  C,  |3 

c  yp 


-mcos  a 


-Cl  (oc+a„,  )vo5>l  +ng^C,„  |3 


c  yp 

Ow 
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cos  (|)v  = 


Cl/^^+Wo,  )'+  ^C„  P 


c  yp ' 

Ow 


T - Vrncos  a  Ci  (a +a„  VoJm+ng — 3 

I  qv  l_ 


Again,  cos  ^  (j)^  +  sin  ^  (j)^  =  1 


Cl  "(a+a„  )'+  3 


ol  '  c  yp  ' 
l  '^w 


mcosa  ((ngfC  "(a+a  )"+  V03,j,^C  B  + 


c  yp ' 

'^W 


s  S  ^  S 

+  2CLja+a„>gM!3^^C,^p+CL/((x+a„jMvci3^P+  ng^C,^  P  -2ClJ(x +a„>gVa3^  p> 

Ow  J 

Simplifying, 


C^^(a+a„J^+  |cc^^P  =  ^mcosa  +  a„J^  +  P  {(n# +(V03,r} 

l^j-w  Jj  j  JJ 

CL/(a+aoJ"+  ;|^Cy^p  =  ^^mcosa  [(ngf  +  (Vo]^  f  ] 

i  j  I 


{  \2  1  (  V- 

la+a„  )  = - -<  — ^mcosa 

^  C,  "  qV^ 


ngf  + (VO] ,,)"]-  |cc^^P 
^  w 


a=—  y^mcosa 
C.,  W  qV^ 


ngf+(Vra,f]-f|^C,  b1  -a„^  (Eqn.l7) 


c  yp  ‘ 

'Jvj 


Substitute  a  into  cos  (|)v 


cos(|)  ^  =  ■ 


1  V 


“  c,  U  ^ 


i  [(ng)"  +  (va3,p  p]-  -^CypP 


-a„  +a„  +  — ^C, 

ol  °l  q  yp 

O  yjj 


mcosa  C 


1  fx 


Cl„  V  ^ 


ng)^+(V03,p)^]-  -^CypB  fig+VtOq, -^CypB 
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which  simplifies  to 


(Eqn.  15) 

Substitute  a  into  sin  (|)v 


118 


(Eqn.  16) 


(Eqn.  18) 


Trim  Condition  9. 

Inner- loop  equations  were  defined  above.  The  parameters  weapon  depression  angle, 
height  above  the  target  and  orbit  radius  were  defined  in  Eigure  10  in  Chapter  3.  The 
equations  describing  the  orbit  parameters  foUow 
V" 

E  =- 
g 

(f)^  =  arccos 
h  =  q  tan((l)^  -i-e) 
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which  combine  to  produce  an  equation  for  the  height. 


f  f  n 

h  = — ,  tan  arccos  —  +e 


gvn  -1 


1  1  1. 

1  -  +  —  tane 


Vn^-1 

1  -  -1  tane  § 


(Eqn.  38) 


To  find  the  maximum  height,  Equation  34  was  differentiated  with  respect  to  n  since  8  was 


constant 


0  =  -(n^  - 1)  '^ntanell- Vn^  -1  tane)+  1+  ^  =tane 


which  was  simphfied. 

n  - 1  n  -1 


n"-l 


Ilhmax 


h_  = 


1  +  tane 
1-tane  g 


(Eqn.  39) 


(f)^  =  arccos  —=  =  45 ' 

*  max  /  /^ 

>  V  Z  , 
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APPENDIX  B:  MATLAB®  ROUTINE 


This  appendix  includes  a  copy  of  the  Matlab®  routine  used  for  the 
implementation  of  the  open- loop  model.  The  control  surface  deflections  and  engine 
settings  were  determined  by  varying  the  input  parameters. 


%  Open-Loop  Control  Laws  for  AFIT 
%  Capt  Gary  Miller 

%  AFIT  2004 

clear  allele 
format  compact 

o. 

o 

%  EQUILIBRIUM  CONSTANTS 
%  Determined  by  specific  flight 
%W=  12712; 

W=  13500; 

S_w=231 . 8; 
b=39.5; 

c_bar=9;  %estimate% 

S_t  =  7  6.5; 
rho0=0 . 0023769; 
g=32 . 174; 

R=1718; 
m=W/ g; 


Thesis 


conditions 

%A/C  weight  at  equilibrium  [lbs] 

%A/C  weight  at  equilibrium  [lbs] 

%A/C  wing  area  [ft^2] 

%Wing  span  of  the  A/C  [ft.] 

%Wing  mean  chord  [ft.] 

%Vertical  Tail  surface  area  [ft^2] 
%S-L  standard  density  [slug/ft^3] 
%Force  of  gravity  [ft/s^2] 

%Universal  gas  const . [ft . lbs/slug*R] 
%Mass  of  A/C  at  equilibrium  [slugs] 


%  OUTER  LOOP 
%  Specify 

%  Beta  or 

%h0=  13879; 
h0=  15000; 
T0=518 . 69-0 . 003 
%V=475; 

%Vo=475; 

V=450; 

Vo=450; 
psi_dot0  =  0 . 0 ; 
phi_dot0  =  0 . 0 ; 
theta_dot0  =  0 . 0 ; 
rho=rho0*exp (-g 
gamma=0 . 0 ; 


CALCULATIONS 
V/Vo,  gamma  and 
phi_v  must  also 
%specif ied% 
%specif ied% 
56616*h0; 

%specif ied% 
%specif ied% 
%specif ied% 
%specif ied% 
%specif ied% 
%specif ied% 
%specif ied% 
*h0/ (R*T0) )  ; 

%specif ied% 


q_bar=rho* (Vo^2) *S_w/2; 
C_y_beta=-0 . 0133*180/pi; 
C_L_alpha=0 .086694*180/pi; 


omegaO  from  flight  conditions, 
be  specified. 

%Equilibrium  Height  [ft]  h0=  13879; 
%Equilibrium  Height  [ft] 

%S-L  standard  temperature  [deg.  R] 
%Wind  axes  velocity  [ft/s] 

%Wind  axes  equil.  velocity  [ft/s] 
%Wind  axes  velocity  [ft/s] 

%Wind  axes  equil.  velocity  [ft/s] 
%Equilibrium  turn  rate  [rad/s] 
%Equil.  roll  angle  change  [rad/s] 
%Equil.  pitch  angle  change  [rad/s] 
%Density,  altitude  corr.  [slug/ft^3] 
%Normallized  flight  path  angle  [rad] 
% (defined  up  as  positive) 
%Dynamic  pressure  [lbs.] 

%Side-Force  Coeff.-  sideslip  [1/rad] 
%Lift  coeff. -angle  of  attack  [1/rad] 
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C_L_o=0 . 010544 ;  %Zero-lift  lift  coefficient  [none] 

C_Do=0 . 02350 ;  %Zero-lift  drag  coefficient  [none] 

C_D_alpha=0 . 00200*180/pi;  %Drag  coeff. -angle  of  attack  [1/rad] 

flag0=l;  %FlagO=l,  beta  specified;  else  phi 

%  Beta  specified  equations 

if  flag0==l;  beta=-0.0;  %Outer  loop  eqns .  w/  beta  specified 

phi=asin (q_bar* ( (V/Vo) ^2) * (C_y_beta*beta*S_t*g/S_w- 
psi_dotO*V*sqrt ( . . . 

( (g^2+ (psi_dotO*V) ^2) * (m* ( (Vo/V) ^2) *cos (gamma) /q_bar) ^2) - . . . 
(C_y_beta*beta*S_t/S_w) ^2) ) / (m*cos (gamma) * (g^2+ (psi_dotO*V) ^2) ) ) ; 

%Roll  angle  point  mass  eqn.  [rad] 

alpha= (sqrt ( (g^2+ (psi_dotO*V) ^2) * ( ( (Vo/V) ^2) *m*cos (gamma) /q_bar) ^2- . . . 
(C_y_beta*S_t*beta/S_w) ^2) -C_L_o) /C_L_alpha; 

%Angle  of  attack  pt .  mass  eqn.  [rad] 
mu=q_bar* ( (V/Vo) ^2) * (C_Do+C_D_alpha* (sqrt ( (g^2+ (psi_dotO*V) ^2) * (m* . . . 

( (Vo/V) ^2) *cos (gamma) /q_bar ) ^2- (C_y_beta*beta*S_t /S_w) ^2 ) - . . . 
C_L_o) /C_L_alpha) /W+sin (gamma) ;  %Engine  throttle  setting 

%  Phi  specified  equations 

else  phi=0.0;  %0uter  loop  eqns.  w/  phi  specified 

beta=- (Vo^2) *W*S_w*cos (gamma) *sin (phi) / (C_y_beta*S_t*q_bar* (V^2) ) ; 

%Sideslip  angle  pt .  mass  eqn.  [rad] 
alpha= ( ( (Vo/V) ^2) *W*cos (gamma) *cos (phi) /q_bar-C_L_o) /C_L_alpha; 

%Angle  of  attack  pt .  mass  eqn.  [rad] 
mu=q_bar* ( (V/Vo) ^2) * (C_Do-C_D_alpha*C_L_o/C_L_alpha) /W+. . . 

C_D_alpha*cos (gamma) *cos (phi) /C_L_alpha+sin (gamma) ; 

%Engine  throttle  setting 

end; 

if  alpha>0.35;  display (' Angle  of  Attack  limit  exceeded');  end; 
if  mu<0; 

display (' Minumum  Throttle  setting  exceeded,  cannot  maintain 

equilibrium' ) ; 

end 

if  mu>l;  display (' Maximum  Throttle  setting  exceeded');  end; 

o. 

o 

%  CONVERT  A/C  wind  axis  angles  to  body  axis  angles 
psi=0.0;  %specified%  %A/C  wind  heading  angle  [rad] 

c_ll=cos (gamma) *cos (psi) *cos (alpha) *cos (beta) + (sin (phi) *sin (gamma) * . . . 
cos (psi) -cos (phi) *sin (psi) ) *cos (alpha) *sin(beta)-(cos(phi)*.  .  . 
sin (gamma) *cos (psi) +sin (phi) *sin (psi) ) *sin  (alpha) ; 

%lst  row,  1st  column  term  matrix 

c_21=cos (gamma) *sin (psi) *cos (alpha) *cos (beta) + (sin (phi) *sin (gamma) * . . . 
sin (psi) +COS (phi) *cos (psi) ) *cos (alpha) *sin(beta)-(cos(phi)*. .  . 
sin (gamma) *sin (psi) -sin (phi) *cos (psi) ) *sin (alpha) ; 

%2nd  row,  1st  column  term  matrix 

c_31  = 

-sin (gamma) *cos (alpha) *cos (beta) +sin (phi) *cos (gamma) *cos (alpha) * . . . 
sin (beta) -cos (phi) *cos (gamma) *sin (alpha) ; 

%3rd  row,  1st  column  term  matrix 
c_32=sin (gamma) *sin (beta) +sin (phi) *cos (gamma) *cos (beta) ; 

%3rd  row,  2nd  column  term  matrix 

c_3  3  = 

-sin (gamma) *sin (alpha) *cos (beta) +sin (phi) *cos (gamma) *sin (alpha) * . . . 
sin (beta) +cos (phi) *cos (gamma) *cos (alpha) ; 

%3rd  row,  3rd  column  term  matrix 
if  abs (-c_31) <pi/2;  theta=asin (-c_31) ; 
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elseif  (abs (sqrt (c_11^2+c_21^2) ) <pi)  & 

(abs (sqrt (c_11^2  +  c_21^2) ) >pi/2)  ; 

theta=pi-acos (sqrt (c_11^2+c_21^2) ) ; %A/C  pitch  angle  <nav.  axes>  [rad] 
end; 

if  abs (c_2 1 / sqrt (c_l 1 ^2+c_2 1^2) ) <pi/2 ; 
psil=asin (c_21/ sqrt (c_11^2+c_21^2) ) ; 

elseif  (abs (c_ll/sqrt (c_11^2+c_21^2) ) <pi)  &  (abs (c_ll/sqrt (c_11^2+ . . . 
c_21^2) ) >pi/2) ; 

psil=pi-acos (c_ll/sqrt (c_11^2+c_21^2) ) ; 
if  psil>pi;  psil=psil-pi;  end;  %Quadrant  check 

end;  %A/C  heading  angle  <nav.  axes>  [rad] 

if  abs (c_32/sqrt (c_11^2+c_21^2) ) <pi/2; 
phil=asin (c_32/ sqrt (c_11^2+c_21^2) ) ; 

elseif  (abs (c_33/sqrt (c_11^2+c_21^2) ) <pi)  &  (abs (c_33/sqrt (c_11^2+ . . . 
c_21^2) ) >pi/2) ; 

phil=pi-acos (c_33/sqrt (c_11^2+c_21^2) ) ; 
if  phil>pi;  phil=phil-pi;  end;  %Quadrant  check 

end;  %A/C  roll  angle  <nav.  axes>  [rad] 

o. 

0 


f lagl=l ; 
if  flagl==l; 

%  CALCULATE  ROTATION  RATES  FROM 
P=-psi_dotO*sin (theta)  ; 
Q=psi_dotO*sin (phil) *cos (theta)  ; 
R=psi_dotO*cos (phil) *cos (theta)  ; 
elseif  flagl==2; 

%  CALCULATE  ROTATION  RATES  FROM 
P  =  0; 

Q=theta_dotO*cos (phil)  ; 
R=-theta_dotO*sin (phil)  ; 
else 

%  CALCULATE  ROTATION  RATES  FROM 
P=phi_dotO ; 

Q=0; 

R=0; 

end; 

%  EQUILIBRIUM  CONSTANTS 
%  AIRCRAFT  MOMENTS  OF  INERTIA 
I_xx=  20000; 

I_yY=  22900; 

I_zz=  40000; 

I_xz=  1980; 


%Define  Rotation  Axis 
%=l>Yaw;  =2>Pitch;  =3>Roll 

EULER  ANGLES  (STEADY  TURN) 
%Equilibrium  Roll  Rate  [rad/s] 
%Equilibrium  Pitch  Rate  [rad/s] 
%Equilibrium  Roll  Rate  [rad/s] 

EULER  ANGLES  (PITCH-OVER) 
%Equilibrium  Roll  Rate  [rad/s] 
%Equilibrium  Pitch  Rate  [rad/s] 
%Equilibrium  Roll  Rate  [rad/s] 

EULER  ANGLES  (STEADY  ROLL) 
%Equilibrium  Roll  Rate  [rad/s] 
%Equilibrium  Pitch  Rate  [rad/s] 
%Equilibrium  Roll  Rate  [rad/s] 


%Mom.  of  Inertia  x-axes  [slug*ft.^2] 
%Mom.  of  Inertia  y-axes  [slug*ft.^2] 
%Mom.  of  Inertia  z-axes  [slug*ft.^2] 
%Prod.  of  Inertia  x/z  [slug*ft.^2] 


%  AIRCRAFT  STABILITY  DERIVATIVES 
C_y_delr=  0 . 002483*180/pi; 

C_y_dela=  -0 . 00107*180/pi; 

C_l_beta=  -0 . 00165*180/pi; 

C_l_dela=  -0 . 00138*180/pi; 

C_l_delr=  0 . 000346*180/pi; 

C_l_p=  -0 . 44634*180/pi; 

C_l_r=  0 . 387836*180/pi; 

C_n_beta=  0 . 001449*180/pi; 

C_n_dela=  -0 . 00002*180/pi; 

C_n_delr=  -0 . 0009*180/pi; 

C_n_p=  -0 . 06415*180/pi; 

C_n_r=  -0 . 11635*180/pi; 


%Side-Force  Coeff, 
%Side-Force  Coeff, 
%Roll  Mom.  Coeff. - 
%Roll  Mom.  Coeff. - 
%Roll  Mom.  Coeff. - 
%Roll  Mom.  Coeff. - 
%Roll  Mom.  Coeff. - 
%Yaw  Mom.  Coeff. - 
%Yaw  Mom.  Coeff. - 
%Yaw  Mom.  Coeff. - 
%Yaw  Mom.  Coeff. - 
%Yaw  Mom.  Coeff. - 


-rud.  def. 
-ail.  def. 

■  sideslip  [ 

■  ail.  def. 

■  rud.  def. 

■  Roll  Rate 

■  Yaw  Rate  [ 
sideslip  [1 
ail.  def.  [ 
rud.  def.  [ 
Roll  Rate  [ 
Yaw  Rate  [s 


[ 1 /rad] 
[ 1 /rad] 
1 /rad] 

[ 1 /rad] 
[ 1 /rad] 
[ s /rad] 
s /rad] 
/rad] 

1 /rad] 

1 /rad] 
s /rad] 
/rad] 
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C_y_p=  0 . 149213*180/pi;  %Side-Force  Coeff.-Roll  Rate  [s/rad] 

C_y_r=  0 . 827684*180/pi;  %Side-Force  Coeff.-Yaw  Rate  [s/rad] 

C_l_pq=  (I_xz* (I_zz+I_xx-I_yy) / ( I_xx* I_z z- I_xz ^2 ) ) ; 

%Roll  Moment  Coeff.  due  Pitch/Roll 
C_l_qr=  ( I_z  z  * ( I_yy-I_z  z ) - I_xx^2 ) / ( I_xx* I_z  z- I_xz  ^2 )  ; 

%Roll  Moment  Coeff.  due  Pitch/Yaw 
C_n_pq=  ( I_xx^2- I_xx* I_yy+ I_xz ^2 ) / ( I_xx* I_z z- I_xz ^2 ) ; 

%Yaw  Moment  Coeff.  due  Pitch/Roll 
C_n_qr=  (I_xz* (I_yy-I_xx-I_zz) ) / ( I_xx* I_z z- I_xz ^2 ) ; 

%Yaw  Moment  Coeff.  due  Pitch/Yaw 


C_mo=  0.052925; 

C_m_dele=  -0 . 01472*180/pi; 

C_m_q=  -15 . 0857*180/pi; 

C_L_q  =  0 . 65200*180/pi; 

C_L_de  =  0 . 006337*180/pi; 

C_m_alpha  =  -0.01557*180/pi; 
C_m_alpha_dot  =  -0 . 09800*180/pi; 
C_d_q  =  0 . 08600*180/pi; 

C_m_p2r2=  I_xz/I_yy; 

C_m_pr=  (I_zz-I_xx) /I_yy; 
c_3=I_zz/ ( I_xx* I_z z- I_xz ^2 ) ; 
c_4=I_xz/ ( I_xx* I_z z-I_xz ^2 ) ; 
c_7=l/I_yy; 

c_8= ( (I_xx-I_yy) * I_xx+ I_xz ^2 ) / (] 

c_9=I_xx/ ( I_xx* I_z z-I_xz ^2 ) ; 

%  THRUST  CHARACTERISTICS 


^Equil.  Pitch  Moment  Coeff.  [none] 
^Pitch  Mom.  Coeff. -elev.  def. [1/rad] 
^Pitch  Mom.  Coeff.  due  Pitch  [s/rad] 
^Lift  Coeff.  due  Pitch[s/rad] 

^Lift  Coeff.  due  elev.  Def.  [1/rad] 
^Pitch  Mom.  Coeff.  due  alpha  [1/rad] 
^Pitch  Mom.  Coeff . -alpha_dot  [s/rad] 
^Drag  Coeff.  due  Pitch  Rate  [s/rad] 
^Pitch  Mom.  Coeff.-Roll  Rate  [s/rad] 
^Pitch  Mom.  Coeff.-  Roll/Yaw  [s/rad] 
^Roll  acc.  eqn.  const.  [ 1 /slug*f t ^2 ] 
^Roll/yaw  eqns .  const.  [ 1 /slug*f t ^2 ] 
^Pitch  acc.  eqn.  const . [ 1 /slug*f t ^2 ] 


\ — 1 

o 

II 

H 

1 

X 

%estimate% 

y_T=3; 

%estimate% 

\ — 1 

O 

II 

H 

1 

N 

%estimate% 

T_max=5000 ; 

%estimate% 

alpha_T=0 . 0 ; 

%estimate% 

T=mu*W ; 

if  T>T_max;  display] 

'Maximum  th. 

per=0 . 5 ; 

%specif ied% 

T_l=per*T; 

T_2=T-T_1 ; 

%  CONTROL  SURFACE 

DEFLECTION  : 

dele_min=  -0.2618; 

dele_max=  0.2618; 

dela_min=  -0.5236; 

dela_max=  0.5236; 

delr_min=  -0.6109; 

delr_max=  0.6109; 

%  TRIM  EQUATIONS 

P_dot=0 . 0 ; 

%specif ied% 

Q_dot=0 . 0 ; 

%specif ied% 

R_dot=0 . 0 ; 

%specif ied% 

alpha_dot=0 . 0 ; 

%specif ied% 

*I_zz-I_ 

X  z  ^  2  )  ; 

%Yaw  acc 

.  eqn.  const. 

[ l/slug*ft 

^2] 

%Yaw  acc 

.  eqn.  const. 

[ l/slug*ft 

-2] 

%x-axis 

distance 

c .  g . 

to  thrust 

[ft 

%y-axis 

distance 

c .  g . 

to  thrust 

[ft 

%z-axis 

distance 

c .  g . 

to  thrust 

[ft 

%Maximum  combined  thrust  [lbs.] 
%Engine  thrust  angle  [rad] 

%Total  engine  trust  [lbs.] 
rust  limit  exceeded');  end; 

%Percentage  of  right  engine  thrust 
%Thrust  from  right  engine  [lbs.] 
%Thrust  from  left  engine  [lbs.] 
LIMITS 

%Negative  Elevator  Def. Limit  [rad] 
%Positive  Elevator  Def.  Limit  [rad] 
%Negative  Aileron  Def.  Limit  [rad] 
%Positive  Aileron  Def.  Limit  [rad] 
%Negative  Rudder  Def.  Limit  [rad] 
%Positive  Rudder  Def.  Limit  [rad] 

%Roll  acceleration  rate  [rad/s^2] 
%Pitch  acceleration  rate  [rad/s^2] 
%Yaw  acceleration  rate  [rad/s^2] 
%Change  in  angle  of  attack  [rad/s] 
%<related  to  theta_dot>% 


c_p=-P_dot+C_l_pq*P*Q+C_l_qr*Q*R+c_3*y_T*sin (alpha_T) * (T_2- 
T_l) +c_4*y_T* . . . 

cos (alpha_T) * (T_2- 

T_l) +q_bar*b* (c_3* (C_l_p*P+C_l_r*R+C_l_beta*beta) + . . . 
c_4 * (C_n_p*P+C_n_r*R+C_n_beta*beta) ) ; 
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%Roll  acceleration  equation  constant 
c_r=-R_dot+C_n_pq*P*Q-C_n_qr*Q*R+c_9*y_T*cos (alpha_T) * (T_2- 
T_l) +c_4*y_T* . . . 

sin (alpha_T) * (T_2- 

T_1 ) +q_bar *b* (c_9* (C_n_p*P+C_n_r*R+C_n_beta*beta) + . . . 
c_4* (C_l_p*P+C_l_r*R+C_l_beta*beta) ) ; 

%Roll  acceleration  equation  constant 
den_lat=q_bar *b* ( (c_3*C_l_dela+c_4*C_n_dela) * (c_9*C_n_delr+c_4 *C_l_delr 

)  .  .  . 

- (c_3*C_l_delr+c_4*C_n_delr) * (c_9*C_n_dela+c_4 *C_l_dela) ) ; 

%Denomenator  lateral  control  eqns . 
dela= (c_r* (c_3*C_l_delr+c_4*C_n_delr) - 
c_p* (c_9*C_n_delr+c_4 *C_l_delr ) ) /. . . 

den_lat;  %Aileron  deflection  equation  [rad] 

delr= (c_p* (c_9*C_n_dela+c_4 *C_l_dela) - 
c_r* (c_3*C_l_dela+c_4*C_n_dela) ) /. . . 

den_lat;  %Rudder  deflection  equation  [rad] 

dele= (c_7* (z_T*cos (alpha_T) +x_T*sin (alpha_T) ) * (T) -C_m_pr*P*R- 
C_m_p2r2  * .  .  . 

(R^2-P^2) -Q_dot) / (C_m_dele*c_7*q_bar*c_bar) - 
(C_mo+C_m_alpha*alpha+ . . . 

C_m_alpha_dot  *alpha_dot+C_m_q*Q) / C_m_dele ; 

%Elevator  deflection  equation  [rad] 
if  dela<dela_min;  display (' Aileron  lower  limit  exceeded');  end; 
if  dela>dela_max;  display (' Aileron  upper  limit  exceeded');  end; 
if  delr<delr_min;  display (' Rudder  lower  limit  exceeded');  end; 
if  delr>delr_max;  display (' Rudder  upper  limit  exceeded');  end; 
if  dele<dele_min;  display (' Elevator  lower  limit  exceeded');  end; 
if  dele>dele_max;  display (' Elevator  upper  limit  exceeded');  end; 
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